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subject:  Page  Revisions,  EoD-TDR-63 -61*5  •  S PrfDC  Program  Document  (AD  1*26  035) 


L  TO=  DDC**  *  a  *  •  < 

I  *  t  * ,  Cameron  Station  *.  *  * 

I  *  »  v'  Alexandria,  Va.  *♦  • 

I  ^ 

i  **  (^  e*  1,  ESD-TDR-63-61*5,  subject,  Special  Perturbation  Weighted  Differen-. 

|  .  tial  Correction  Prefer  am,  was  transmitted  to  DDQ.  on  31  December  1963. 

|  *  •  •  \  •  4  •  . 
i  *  \  o  • 

^  J  2.  The  attached  page  revisions  are  forwarded  for  inclusion  in  the. 
basic  report..  • 
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24  .’une  1964 


SUMMARY  OF  CURREKT  SPWDC  IK PUT  CARD  'FORMATS  .  #  , 

•  .  The  programs  SPWDC  end  dSSWGTE  have  undergone  changes  by  agenc/ies 
•  •  i 

both  at* Colorado  Springs  and  at  Bedford  which  Invalidate  some  of  tne 

input  card  formats  as  described, in  the  program  document.  Althqjigh  these 

changes  have  been  'documented  through  program  change  notes,  some  of  the  *. 

program  users  have  not  become  aware  of  the* changes  in  prdgram  operation 
*  ,  •  «  # 
because  of  the  limited  distribution  ofchangq  notes.  In  order  to  sum* 

marize  .the. format  changes  in  a* single  source  and  to  provide  for  ade¬ 
quate  distribution,  this*  m ea-j  has  fcejn  prepared.  Figures  shqwlng  'the 
revised  formats  for  the  changed  cards  are  included.  It, is  suggested 
that  these  be  inserted  into  each  copy  of  the  s'PWDC  Program  Document,* 
(ESD-TDR-63-645) .  *  *  .*  * 


*  With  respect  to  observation  precession,  It  should  be*noted  that 

if  column  80  of  an  observation  card  is  blank,  and  it  is  an  optical ’type,, 

(right-ascension,  declination)  then  it  is  precessed  in  accordance  with. 

^  •  •  •  •  • 
the  standard  system  technique.  'That  is  to  say,  that  if* right-ascension 
•  •  *  » 
and  declination  are  observed,’  and  the  equipment  type  is  not  16,'  then  it* 

••  •  ,  .  •  •  •  . 
is  precessed  from*  the  year  which  Is  indicated  in,  column  70.  co  the  present 

and  if  it  is  a  field  reduced  Baker-Kunn  observation  (equipment  typo  16), 

«  ,  •  •  *•  . 

•it  is  precessed,  from  1855  to  the  present  if  the  declination  is  greater  .* 

than  -22°  ani  from’f875  if  ft  is  less.  ’  •  ** 


*  Prepared  for  the  496L  SPO  by  Aeronutronic, 
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Reference:  CH  X-18B  with  sttacliaents,  CN  T-335 
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*  •  •  •  *  « .  •  •  - 

^  ••  «*•*  •*  •  • 

|  •  Field,  •  Colunn  .  ^onte^ts  *  a  a  *  a  •  Code  .  , 

*  *  *  *  *.  ,  a  ,  a*  ,  ,  *  .  »  .  .  ..  .  .  •  » 

•  *  *  *  i  .  .  1  -  8  a  ."SPWDC&AA"  •  *  *  . 

#  .  •  ••  •  *  *« 
I  •  v2  *  9a-#ll  Satellite  Number  (optional)  •  •  *  '  . 

c  *-  #  •  •  •  •  •  •  *  •  #»*'*,  * '  ** 

_  *  3  t  12  Auxiliary  veight  tape  input* (2«no,  0»yes)  ’  •*  • 

^  •  •  t  *  |  ^  ^  j  *•«»*•** 

•  3*  •  t>  *  *  «.  *13  ,  *  P.iw  data  observation  tape  Input  (0»no,  f-yoe)  .  * 

*1*  ,  *  a  •  *  *  •  *  •  • 

*  5  .  14  -Mim  fJalunvA  be  blank  for  SPUDC.  •  • 


5  ,14  Jlusti  always  be  blank  for  SPWDC.  ••  ’  •  .  *  *  • 

6  *  15  *  Print  Input  observations  prior  to  lifc.  (0»no*,  l-yes)* .  .  •  • 

•  •  }  ae»  *  •  .  . 

7  16  -  76  *  Not  used.  »  •  4 

a  *  *  a  •  .  .  .  ® 


o 

I  .  • 

,  »  ••  10* 


77  Weighting  flag.  (l*mo*,  0*yes)  *  * 

78  -  79,  Always  "A0’\  >  *  *(  • 

80  Always  "P"  **  *  •  * 

•  ,• a  o  * 
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FIGURE  17.  0BS&GTE  PARA1ETER  CARi) 
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1  Field  * Column  *  Content#  ’ 

V  •  6  ' 

'«  1  1  tit  Integration  made  (  A"vartable,  1-fixed) 

2  2  -  11  '  .  At  (ftlnutes)  (positive  makes  correction  epflcb 

.  *  -.before  obs.,  -negative  makes  It  after'obs,) 

*  ♦  *  • 

. .  3  .  »*.  12  Bulge  perturbation  flag 

’4  *  •  !  13  Drag  oerturbatlon  flag 

:  ...  » 

£»  .*14.  Radiation  pressure  perturbation  flag 

•  * 

" '.6‘'  ’  *15  ’  New  epoeh  mode  (0*rev*  l»time*,*  2»last  ob.) 

«  ’  *  *  m  *  ^ 

7  *J16  -  29  •  New  epoch  (day  number  "in  year  and  fraction*  or 

«  resolution  number.)* 

•*  »  «  a  ••• 

r‘fi  :  30-  -  36  Elements  to  correct:  n,  a^,  ayN,  0^,  XJ,  1,  m» 

9  37 .»  38  *  Number  of  times  to  repeat  correction. > 

#  •  ••  *  *  * 

•  \0  39  Aq  chick  flag  (used  only  if  or  7  ^foments  are 

•  ®  .  *  *  ’’  being  corrected.) 

•  t  •  • 

•  11  *  40c-  47  Maximum  &q  (km)  (used  if  col.  39  1.) 

#12  *  ;  48  -  55  #  ABSMX  (km)  #  # 

•  13  *•  36  -  6*3  *  ABSMX 2*  (ka/sec)  .  . 

*/  .  •  •  .  .  . 

14.  £4  -  71*  n  (rms  multiplier  for  rejection.) 


*  FIGURE  19.  DIFFERENTIAL  CORRECTION  CONTROL  CARD,  Part  1  of  2 
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Concents  • 

Code 

15 

•  ^  72 

Residugl  output  flag  (0-never,  l»first  and  last 

• 

passes  only,  2»all  passes.) 

2 

16 

73 

,  Specif  "n  and  U"  correction  on  first  pass. 

4 

\1 

?4  -  77 

Relative  change  in  successive  rns's  to  determine 

• 

convergence,  (if  blank,  .05  is  used.) 

3  • 

18 

78  -  79 

Always  "tu'l 

• 

19 

9 

*  30 

9 

Always  ''P”.  * 

•  * 

* 

• 

**  * 

FIGURE  19.  DIFFERENTIAL  CORRECTION  CONTROL  CARD,  P*rt  2  of  2  • 
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Field 

1 

2 

3 

4 

• 
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7 

8 
9 

10 

11 

12 


Colusa  Contents  Code 

1  nt. integration  nwA  (  A -variable,  Infixed)  7 

o  |  9 

2-11  At  (t$lnutes>  (sign  determined  from  7?  card.)  3 

12  Bulge  perturbation  flag  f>  4 

* 

13  Brag  perturbation  flag  *  4 

a 

9  14  Radiation  pressure  perturbation  flag  *  4 

e  • 

,15  -  16  Print  rflag  j  1  -  t,  r,  i  *\  Add  flags  for 

»  I  2  -  t,  a,  e,  1,  0.  u),  U  V  combination 

*  14  -  t,  0o  Ag,  h  of  of  data. 

1  W  Binary  ephemrtis  ta^  ovt-nt  flag 

18  Prediction  reliability  output  flltg 

• 

19  Punched  cards  (t,  r)  (Also  leaves  data  In  core 

for  later  use  by  GIP^R)  a  •  • 

ft 

20  ••  7^  Not  used.  * 

73-  79  9  Always  "16".  *  *  *  # 

80  Alv«y«#,,P,,a 


FIGURE  22.  TIKE  PREDICTION  CARD  ti 
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Assumed  Decimal  Points 


mill 
1  3  3  3  3  3 
4  4  (444 
S»SSSJ 
6  S  tit  I  ft 


IIIIUII 


Column 


Contents 

•  * 

Sensor  number 


1-4  Sensor  number  ®  1 

•  5  -  11  Station  latitude  (deg.  north)  (assumed? dec.  pnt.)  3 

•  *  • 

12  -  19  Station  longitude  (deg.  west)  (asserted  del.  pnt.)  3 


12  -  19 
20  -  25 
26  -  77 
78  -  79 


Station  altitude  (meters)  (assumed  dec.  pnt.) 
Not  used. 

Always  "A9". 

Always  "P".  9  • 


Note:  Trailitig»zeros  must  be  punehed^ln  fields  2  and  3. 


-»  FIGURE  25.  STATION  PASS  PREDICTION  CARD  #2 


Reference:  CN*s  T-227,  T-248,  T-276,  T-286,  T-299,  &  T-313. 
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Column 


4-77 
78  -  79 


Content*  _  •  Coiu 

€>  • 

• 

DC  Flag  #4  • 

Prediction  (0-no) (1-Tlme) (2~Statlon)  2 

0#or  blank  -  Print  only  (at,  U,  end  /J  for  "n  and  2 
U"  correction  pas*. 

1  -  Print  also  Ap,  />cosSA<*,phS  ,  Ap  for  "n  and 
U"  correction  pasa. 

2  ’  Print  Ap,  cosSAtf,  AS,  A/5 ,  In  degree*  for 

all  phases.  *  • 

3  -  Hake  on one  pats.  Additional  realduals  In 

U,V,W  coordinates  are  computed  and  printed,  a* 
veil  as  residual  means  and  standard  deviations 
for  each  of  the  sensors.  * 

4  -  Like  option  3  except  that  residuals  In  oblerved 

quantities  are  punched  on  cards  as*uell  as  printed. 
(Format  given  in  CN  T-313.)  • 

•Not  used, ■ 

Always  "10". 

Always  "P" .  • 


FIGURE  26.  PROGRAM  EXECUTION  CARD 
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Reference:  CH  R-188  with  attachments. 
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•  Field 


Column 


Contents 


•  1 

1-3 

Sensor  number 

1 

2 

4-5 

Always  "11", 

J 

«  6-8 

Not  used.  . 

4 

9-16 

Always  "CSGBIASA"  # 

5 

17  -  24 

Range  standard  deviation  (5cm)  •  * 

3 

6 

a  25  -  32  a 

Range  bias  (km)  .  s 

3 

7 

33  -  40 

Range-rate  standard  deviation  (km/sec) 

3 

8 

41-48 

Range-rate® Vas  (1cm/  sec)  *  • 

3 

9 

49-56 

Aslmuth  (or  rt.  asc.)  standard  deviation  (deg) 

3 

Q  10 

57  -  64, 

Aslmuth  (pr  ft.  asc.)  bias  (deg) 

3 

11 

65  -  72 

Elevation  (or  decl.)  standard  deviation  (deg) 

3 

12 

73 

Multiple  punch  (11-8-2) 

13 

74  -  80 

• 

Not  used,  # 

*  FIGURE  30A 

.  CSGBIAS  FUNCTION  WEIGHTING  CARD  #1* 

i  « 

• 

• 
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'.Reference:  CN  K-188  with  attachments. 
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Column  Contents  9 

1-3  Sensor  number 

4-5  Always  "10". 

6-8  Not  usec<. 

9-16  Always  "CSGBIASA" 

17  -  24  Elevation  (or  decl.)  bias  (deg) 

25  -  32  Observation  time  bias  (sec) 

33  -  40  Not  used 

41  Multiple  punch  (11-8-2) 

42  -  80  Not  used 


FIGURE  30B.  CS^BIAS  FUNCTION  WEIGHTING  CARD  #2 
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'•  Reference:  CN  K-112  with  attachment. 
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Cod*-, 


i 

1-3 

2 

4-5 

3 

6-8 

4 

9-16 

5 

17  -  24 

6 

25 

-  32 

7 

33 

-  40 

8 

41 

-  48 

9 

49 

-  56 

10 

57 

-  64 

11 

65 

-  72 

12 

73 

-  79 

13 

80 

Sensor  number  1 

Always  "10".  •  # 

Hot  used.  , 

Always  "SIGTABU"  •  • 

Alternate  value  selector  (anywhere  in  field)  2 

0  or  A  -  Primary  value 

1  -  1st  alternate  value  if  compiled  in  OBSWGTE 

2  -  2nd  alternate  value  If  compiled  in  OBSKGTE 
Scale  factor.  All  tabular  values  of  itandard  deviations  3 
are  multiplied  by  this  factor.  (If  blank,  a  value  of 


one  is  assumed.) 

Range  bias  (km)  9  3 

Elevation  (or  decl.)  bias,  (deg)  3 

Azimuth  (or  rt.  esc.)  bias,  (deg)  3 

Range-rate  bias,  (km/sec)  3 

Observation  time  bias,  (min.)  3 

Not  used. 

Multiple  punch  (11-8-2) 


FIGURE  30C.  SIGTAB1  FUNCTION  V'EIGKtINC  CARD 
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Whep  US  Government  drawings,  specifications  or  other  data  are 
used  for  any  purpose  other  than  a  definitely  related  government 
procurement  operation,  the  government  thereby  incurs  no  respons¬ 
ibility  nor  any  obligation  whatsoever;  and  the  fact  that  the 
government  may  have  formulated,  furnished,  or  in  any  way  supplied 
the  said  drawings,  specifications,  or  other  data  is  not  to  be 
regarded  by  implication  or  otherwise,  as  in  any  manner  licensing 
the  holder  or  any  other  person  or  corporation,  or  conveying  any 
rights  or  permission  to  manufacture,  use,  or  sell  any  patented 
invention  that  may  in  any  way  be  related  thereto. 
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ABSTRACT 


This  report  describes  a  computer  program  which  differentially 
corrects  and/or  predicts  the  orbit  of  a  geocentric  satellite.  The  dynamics 
of  the  program  are  based  on  a  variation  of  parameters  formulation  with  the 
perturbative  acceleration  being  numerically  integrated.  A  weighted  correc¬ 
tion  process  is  incorporated  which  considers  the  relative  accuracy  of  the 
observational  data  in  a  least-squares  fit.  Included  in  the  document  are 
the  program  description,  complete  formulation,  program  operating  instructions, 
flow  diagrams,  and  test  cases. 
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SECTION  1 


INTRODUCTION 


Aeronutronic  has  developed  an  operational  weighted  differential 
correction  program,  using  the  special  perturbations  variation  of  parameters 
approach.  Operating  under  the  496L  B-2  Semi-automatic  Programming  System, 
the  program  will  accept  sensor  observations,  assign  weights  to  the  data, 
perform  a  differential  correction  and  predict  future  position  and  velocity. 

The  purpose  of  the  Special  Perturbations  Weighted  Differential 
Correction  (SPWDC)  Program  is  to  provide  highly  precise  predictions  of 
satellite  position.  The  operational  capability  to  utilize  weights  when 
performing  a  differential  correction  in  addition  to  the  highly  accurate 
Special  Perturbations  formulation  will  provide  far  more  accurate  position 
prediction  than  the  Simplified  General  Perturbations  approach.  The  latter 
concept  was  designed  to  provide  sufficient  accuracy  to  maintain  surveillance 
on  all  earth  satellites  with  minimum  computer  time  requirements.  This  sur¬ 
veillance  mission  requires  only  such  accuracy  as  "o  ensure  acquisition  and 
recognition  of  every  satellite  by  sensors  capable  of  observing  it. 

• 

The  SPWDC  Program,  therefore,  supplies  a  means  for  correcting 
orbital  parameters  and  predicting  in  an  extremely  accurate  fashion  for 
missions  requiring  precise  position  determination.  Operationally,  it 
functions  in  the  B-2  System  running  in  the  Schedule  Tape  Mode. 


Section  2  of  this  document  is  a  description  of  the  SPWDC  Program. 
The  formulation  by  which  the  program  generates  the  numerically  integrated 
position  is  given  in  Section  3  of  the  report.  Operating  instructions  are 
detailed  in  Section  4,  with  parameter  card  formats  shown  in  Appendix  I.  The 
main  flow  of  the  program  is  diagramed  in  Section  5,  and  detailed  flow  data 
in  Appendix  E.  Compiled  test  cases  are  outlined  in  Section  6.  A  glossary 
of  terms  will  be  found  in  Appendix  III. 


SECTION  2 


PROGRAM  DESCRIPTION 


The  SFWDC  Program  contains  three  main  functions.  These 
functions,  shown  schematically  in  Figure  1,  are:  OBSEPvVATION  WEIGHTING, 
DIFFERENTIAL  CORRECTION,  and  PREDICTION.  Each  of  these  functions  is 
iescribed  here,  in  terms  of  its  general  features,  required  parameters, 
and  function  controls  both  internal  and  output.  Required  parameters 
are  those  numerical  quantities  describing  the  physical  conditions  to  be 
simulated  in  the  program.  Internal  controls  are  those  quantities  which 
specify  the  program  logic  to  be  executed,  and  output  controls  are  those 
quantities  prescribing  the  choice  of  available  output. 


2.1  OBSERVATION  WEIGHTING 

The  Observation  Weighting  function  of  the  program  is  to 
reformat  input  observations,  assigning  to  each  observed  quantity  its 
standard  deviation.  These  accuracy  estimates  are  provided  to  the 
routine  in  three  ways: 

(1)  By  assigning  a  constant  accuracy  estimate  for  all 
observations  from  a  given  sensor; 

(2)  By  assigning  an  accuracy  estimate  to  each  observa¬ 
tion  being  used;  or 

(3)  By  indicating  a  programmed  weight  function  which  is 
used  to  compute  the  observation  parameters. 

The  assignment  of  weight  factors  to  the  observations  may  be  made  by  any 
combination  of  the  above  methods,  with  the  restriction  that  any  given 
sensor  will  use  only  one  of  three  choices. 
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FIGURE  I.  SPWDC  PROGRAM  FUNCTIONS 


After  the  complete  set  of  observations  has  been  examined  and 
each  weight  assignment  has  been  made,  the  observations  are  ordered  chrono¬ 
logically.  This  feature  is  necessary  in  the  Observation  Weighting  function 
in  order  to  ensure  an  efficient  correction  process.  Since  the  correction 
program  must  represent  the  observations  at  each  actual  observation  time  by 
numerically  integrating  to  that  time,  it  is  highly  desirable  that  the 
observations  be  ordered  in  time  so  the  integration  takes  place  entirely  in 
one  direction. 


a.  Observation  Weighting  Parameters 

Observations  are  supplied  by: 

(1)  Standard  SPS  Observation  Card,  or 

(2)  SRADU  tape. 

Weighting  information  is  supplied  by  the  Weighting  Tape  (see 

Section  4) . 

b.  Observation  Weighting  Controls 

Internal  controls  are  as  follows: 

(1)  Weighting  information  may  be  picked  up  either  from  the 
Weighting  Tape,  or  a  set  of  accuracy  estimates  may  be 
compiled  into  the  Observation  Weighting  Program  in  the 
weight  by  sensor  mode.  Normally,  the  former  selection 
will  be  made;  however,  it  may  be  to  the  user's  advantage 
to  program  the  sensor  weight  values  for  automatic  inter¬ 
nal  assignment. 

(2)  Weighting  data  may  or  may  not  be  assigned  to  the 
observations.  It  should  be  noted  that  if  weights  are 
to  be  assigned,  they  must  be  assigned  to  all  observa¬ 
tions.  It  is  not  possible  to  mix  weighted  and  unweighted 
data . 

(3)  A  selection  of  weighting  mode  is  required  to  indicate 
the  choice  of  weighting  by  sensor,  by  observation,  or 
by  function. 


-5- 


No  output  controls  are  available  to  the  Observation  Weighting 


function. 


2.2  DIFFERENTIAL  CORRECTION 

The  Differential  Correction  function  accepts  the  reformatted 
observations  from  the  Observation  Weighting  function  and  performs  a 
weighted  or  unweighted  differential  correction  on  the  orbit  element  set. 
Ephemeris  computation  is  carried  out  in  a  special  perturbations  variation 
of  parameters  formulation.  This  formulation  numerically  integrates  the 
perturbative  accelerations  influencing  the  parameters  of  a;,  instantaneous 
two-body  reference  orbit. 

The  numerical  integration  is  performed  by  a  fourth-order  Runge- 
Kutta  integration  routine.  A  fixed  step  size  may  be  specified  for  this 
numerical  integration,  or  it  may  be  used  in  a  variable  mode.  In  the 
variable  mode,  the  initial  step  size  is  specified  and  subsequently 
increased  or  decreased  internally,  based  on  error  control  parameters. 

Perturbations  handled  by  the  program  include:  zonal  harmonic 
bulge,  atmospheric  drag,  and  solar  radiation  pressure.  These  effects  may 
be  controlled  to  incorporate  any  combination  in  a  particular  case. 

o 

A  correction  can  be  obtained  for  any  of  the  six  orbit  elements, 
and  may  also  inc!  ide  the  satellite  mass.  This  latter  parameter  actually 
represents  a  correction  on  the  ballistic  parameter  CqA  >  where  is  the 

m 

drag  coefficient,  A  is  the  satellite  cross-sectional  area,  and  m  is  the 
satellite  mass. 

First-order  partials  relating  variations  in  observed  quantities 
to  variations  in  orbit  parameters  are  the  basis  of  the  correction  equa¬ 
tions.  The  correction  for  mass,  or  ballistic  parameter,  uses  an  equation 
derived  on  the  assumption  of  near  circular  eccentricity.  Therefore,  for 
the  general  case,  a  variant  orbit  approach  is  used,  whereby  the  satellite 
mass  is  varied  and  a  complete  integration  is  performed.  This  ephemeris 
is  then  used  to  form  numerical  partials  between  the  observed  quantities 
and  the  variation  in  mass. 

Special  correction  passes  may  be  taken.  For  instance,  it  may 
be  desirable  to  correct  the  mean  motion  initially  before  attempting  a 
complete  correction.  This  will  reduce  the  timing  residual  and  provide 
a  better  fit  on  subsequent  correction  passes.  Also,  an  incernal  check 
has  been  built  in  which  examines  the  change  in  the  perigee  distance 
between  successive  corrections.  If  the  perigee  distance  varies  too 


greatly,  the  correction  is  automatically  repeated  without  correcting 
the  orbital  eccentricity.  This  prevents  a  possible  divergent  correction 
which  may  occur  due  to  non-linear  drag  effects. 

The  correction  process  will  repeat  up  to  a  maximum  number  of 
iterations,  or  until  the  root-mean-square  (RMS)  of  the  accepted  residuals 
has  converged  within  5%  comparing  successive  values.  Accepted  residuals 
are  those  that  are  not  rejected  on  the  basis  of  an  absolute  maximum  or 
a  relative  check  based  on  a  limit  of  1.5  times  the  residual  RMS.  In  a 
weighted  differential  correction  this  is  the  weighted  RMS;  however, 
the  convergence  check  is  still  against  the  unweighted  value.  In 
addition,  weighted  runs  always  perform  an  unweighted  correction  on  the 
first  iteration  in  order  to  fit  the  elements  as  well  as  possible  before 
applying  weight  data.  An  option  is  also  provided,  enabling  the  post 
correction  epoch  to  be  updated  to  a  specific  revolution  number  or  time. 

a.  Differential  Correction  Parameters 


Orbit  elements  are  supplied  by: 

(1)  The  standard  card  format,  or 

(2)  The  SEAI  file. 

Vehicle  data  are  mass,  diameter,  and  reflectivity, 
b.  Differential  Correction  Contro.s 


Internal  controls  are  as  follows: 

(1)  The  integration  mode,  variable  or  fixed,  is 
specified. 

(2)  The  integration  step  size  is  specified.  If  the 
mode  is  variable,  then  this  denotes  the  initial 
step  size. 

(3)  Selection  of  perturbations  to  be  used:  bulge, 
drag,  radiation  pressure. 

(4)  Specification  of  the  new  epoch: 

(a)  Revolution  number  (actual  number,  not  refer¬ 
enced  to  current  epoch) 

(b)  Time  (absolute) ,  or 


(c)  Time  of  last  observation. 


(5)  Specification  of  the  elements  to  be  corrected  Any 
combination  may  be  selected;  however,  in  some  instances 
the  corrections  are  not  independent  and  will  influence 
the  other  parameters. 

(6)  Special  correction  passes: 

(a)  Correct  mean  motion  only  on  first  pass 

•  (b)  Maximum  value  allowed  for  change  in  perigee 

distance 

(7)  Specification  of  rejection  criteria: 

(a)  Absolute  maximum  for  displacement  residuals, 

(b)  Absolute  maximum  for  range-rate  residuals,  and 

(c)  Multiplying  factor  applied  to  the  residual  RMS 
for  relative  rejection. 

(8)  Error  co-  ols  for  Runge-Kutta  integration: 

(a)  Absolute  controls 

(b)  Relative  controls 

One  output  control  is  available  to  the  Differential  Correction 
function.  This  establishes  the  frequency  of  output  of  the  residuals  as: 
never,  first  and  last  set,  or  every  set. 


2 . 3  PREDICTION 

The  Prediction  function  of  the  program  may  be  used  to  obtain  future 
position  and  velocity  data  from  an  input  element  set  or  a  corrected  set  from 
the  Differential  Correction  function.  Two  methods  are  provided  for  the  speci¬ 
fication  of  the  prediction  time.  A  series  of  prediction  points  may  be 
specified  by  time  or  by  the  closest  point  of  approach  to  some  ’ocation  on  the 
earth . 
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Ephemeris  computation  in  the  Prediction  function  is  carried  out 
in  the  same  manner  as  in  the  Differential  Correction,  i.e.,  using  a 
special  perturbations  variation  of  parameters  formulation  where  the 
integration  is  according  to  a  fourth-order  Runge-Kutta  process .  The  mode 
of  integration,  integration  step  size,  and  the  perturbations  to  be  con¬ 
sidered  must  be  supplied. 

In  addition  to  the  prediction,  this  function  has  the  capability 
of  estimating  the  reliability  of  the  prediction  point.  The  reliability 
check  computes  the  variance  in  both  inertial  and  radial,  transverse  and 
orthogonal  components  of  position  and  velocity.  It  is  based  on  the 
variance-covariance  matrix  of  orbit  elements  derived  from  the  Differ¬ 
ential  Correction  function. 

a.  Prediction  Parameters 


Prediction  by  time  requires  the  following  parameters: 

(1)  Initial  prediction  time  (absolute) 

(2)  Prediction  interval  (minutes) 

(3)  Number  of  points  to  be  obtained 

Prediction  by  station  pass  requires  the  following  parameters: 

(1)  A  station  number  identification.  This  would 
pertain  to  the  SEAl  File  designation  if  the 
location  is  to  be  one  of  the  SPADATS  sensors. 

Otherwise,  any  number  designation  will  suffice. 

(2)  Latitude,  longitude  (west),  and  altitude  of 
the  station,  required  if  the  SEAI  File  is  not 
being  used. 

(3)  The  length  of  time  to  simulate  the  ephemeris 

on  either  side  of  the  closest  point  of  approach. 

(4)  The  maximum  elevation  angle  below  which  no 
ephemeris  data  will  be  obtained.  This  takes 
precedence  over  item  (3)  above. 

(5)  The  interval  of  time  to  be  used  in  obtaining 
the  ephemeris . 

(6)  The  number  of  passes  to  be  considered. 
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b.  Prediction  Controls 


Internal  controls  a^e  as  follows: 

(1)  The  integration  mode,  variable  or  fixed,  is 
specified . 

(2)  The  integration  step  size  is  specified.  As 
in  the  Differential  Correction,  if  the  mode 
is  variable,  then  this  denotes  the  initial 
step  size . 

(3)  Selection  of  perturbations  to  be  used:  bulge, 
drag,  radiation  pressure. 

Output  controls  are  as  follows: 

(1)  Print  any  combination  of 

(a)  time,  position,  and  velocity  (t,  x,  y, 
z,  x,  y,  z) 

(b)  time,  osculating  elements  (t,  a,  e,  i, 

0,  ,  CO,  L) 

(c)  time,  ground  track  data  (t,  satellite  lati¬ 
tude,  satellite  east  longitude,  and  satellite 
altitude) 

(d)  time,  look  angles  (t,  range,  range-rate, 
azimuth, elevation)  available  only  for 
station  pass  prediction 

(2)  Print  reliability  estimate  consisting  of  the 
standard  deviation  in  the  position  and  velocity 
predictions  in  inertial  coordinates  and  orbital 
coordinates . 

(3)  Prepare  a  binary  ephemeris  tape  containing  t,  x, 
y,  z,  x,  y,  z. 

(4)  Punch  position  cards  containing  t,  x,  y,  z. 


-10- 


SECTION  3 


FORMULATION 


The  Special  Perturbation  Weighted  Differential  Correction  (SPWDC) 
Program  employs  a  variation  of  parameters  formulation  to  generate  numerically 
integrated  positions.  This  ephemeris  is  computed  using  the  N,  M  parameters, 
namely,  mean  motion,  n  ,  mean  longitude,  L  ;  two  parameters  embodying  eccen¬ 
tricity,  e,  and  argument  of  perigee,  =  e  cos  (D  and  a^No  =  e  sin  co; 

and  a  vector  defining  the  angular  momentum,  hQ  whose  components  are 


hK  =  v/pT 

sin 

sin  i0 

o 

II 

o 

>> 

cos 

sin  *0 

hz  =  n/pT 

cos 

io 

o 


where  p  is  the  semi-latus  rectum,  Q,q  is  the  right  ascension  of  the  ascending 
node,  and  iQ  is  the  orbital  inclination  (see  Figure  2) ,  The  subscript  identi¬ 
fies  these  quantities  as  the  orbital  parameters  at  epoch,  t  .  These  parameters 
do  not  exhibit  singularities  at  zero  eccentricity .  A  parameter  adoption  in 
the  differential  correction  section  of  the  program  effectively  limits  the 
program  application  to  orbits  having  inclinations  greater  than  two  degrees. 
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NORTH 

CELESTIAL  POLE 


FIGURE  2.  PROJECTION  OF  ORBIT  ON  CELESTIAL  SPHERE,  WITH 
ORIENTATION  UNIT  VECTORS  AND  ANGLES  DISPLAYED 


The  formulation  accounts  for  perturbations  due  to  drag  forces  and 
the  primary  effects  of  the  Earth'o  asphericity^*  and  to  radiation  pressure ^  . 
The  elements  LQ,  nQ,  ax^Q,  and  Jiq  vary  from  point  to  point  on  the  resulting 
ephemeris  in  such  a  manner  that  they  satisfy  a  syscem  of  differential  equations 
in  the  dependent  variables,  the  elements,  and  the  independent  variable,  time,  t. 
The  solution  to  the  system  of  differential  equations  is  achieved  in  a  numerical 
integration  process,  executed  by  means  of  the  Runge-Kutta  procedure.  The  pro¬ 
gram  proceeds  to  evaluate  the  elements  from  point  to  point  by  using  either  a 
fixed  grid  size  for  the  variable  t  or  a  variable  grid  size  controlled  by  the 
program. 


The  program  will  function  in  a  differential  correction  mode  and/or 
an  ephemeris  mode.  In  the  differential  correction  mode,  a  given  set  of  elements 
is  improved  and  updated  by  means  of  a  least-squares  fit  to  a  set  of  observations. 
These  observations  may  be  a  weighted  set  or  an  unweighted  set  with  the  correc¬ 
tion  functioning  accordingly.  In  the  prediction  mode,  future  position  and 
velocity  data  are  provided  by  integrating  to  the  specified  point,  using  the 
corrected  set  of  elements  or  an  element  set  introduced  for  the.  purpose  of  pre¬ 
diction  only. 

3.1  OBSERVATION  WEIGHTING 

One  of  the  basic  considerations  to  be  included  in  the  development  of 
accurate  Special  Perturbations  programs  is  that  of  assignment  of  weighting 
factors  to  all  observations  for  subsequent  use  in  the  differential  correction 
process.  Because  of  the  restrictive  core  allocation  imposed  by  the  operational 
system,  the  weight  assignment  is  performed  independent  of  the  remainder  of  the 
processing. 

The  two  main  design  features  of  the  resulting  proprocessing  program 
are  (1)  to  reformat  the  observational  data,  adding  weights  for  all  the  quantities, 
and  (2)  to  order  the  observations  in  time,  for  either  a  forward  or  backward 
integration. 

These  design  criteria  were  incorporated  into  the  observation  weighting 
function  (OBSWGT)  of  SPWDC.  The  functional  flow  between  OBSWGT  and  SPWDC  is 
shown  in  Figure  3.  The  general  specifications  of  the  program  are: 


*Superscripts  indicate  References  at  end  of  document. 
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FIGURE  3.  OBSWGT  FUNCTIONAL  FLOW 


ASSIGN 
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OBSERVATION 


CALL  IN  SPWDC 
AND  TRANSFER 
CONTROL 
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(1)  OBSWGT  will  operate  before  the  Special  Perturbation 
Program.  OBSWGT  will  apply  weights  to  all  the  obser¬ 
vations,  leave  the  observations  in  core,  and  call 
(and  exit  to)  SPWDC  from  the  master  program  tape. 

(2)  OBSWGT  will  operate  as  a  schedule  tape  job.  Standard 
input  will  consist  of  combinations  of  parameter,  element, 
sensor  and  observations  cards. 

(3)  Additional  input  data  will  be  contained  on  an  auxiliary 
tape.  Data  on  this  tape  will  contain  weight  factors  for 
each  sensor  and/or  observation  or  will  contain  controls 
for  determining  the  method  to  be  used  for  computing  the 
weights.  Operationally,  this  will  be  !;nown  as  the 
Weighting  Tape. 

(a)  A  single  set  of  weights  (at  most,  four 
quantities)  can  be  specified  for  a  sensor.  These 
weights  then  apply  to  all  observations  reported 
by  the  sensor. 

(b)  A  single  set  of  weights  can  be  specified  for  each 
observation. 

(c)  A  routine  can  be  specified  by  sensor  to  compute 
weight  factors  for  each  observation  reported  by  the 
sensor . 

(4)  A  subroutine  will  be  made  available  to  the  Special  Pertur¬ 
bation  program  to  retrieve  the  observation  quantities,  weight 
factors,  and  to  retrieve  the  required  sensor  data. 

(5)  The  routines  for  computing  weight  factors  (mentioned  in  (3) 

(c)  above)  will  have  a  standard  calling  format  and  be  con¬ 
tained  as  a  subroutine  in  the  OBSWGT  program.  These  routines 
can  only  be  added,  deleted,  or  changed  by  recompilations. 
However,  the  parameters  used  by  the  routines  can  be  changed 
(specified  as  input  quantities)  prior  to  each  run. 

(6)  All  weighting  procedures  will  be  general  so  that  the  specific 
function  to  be  applied  is  determined  on  the  Weighting  Tape  by 
pairing  a  sensor  with  a  function  as  an  input  quantity.  In 
this  way  routines  can  apply  to  more  than  one  sensor,  reducing 
duplication  of  routines. 

(7)  A  special  bypass  mode  will  be  available  to  process  observations 
without  applying  weight  factors. 
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3.2  VARIATION  OF  PARAMETERS  EPHErfERIS  COMPUTATION 


The  following  paragraphs  detail  the  formulation  of  the  represen¬ 
tation,  accelerations,  and  integration  of  the  orbit.  The  representation 
calculates  position  and  velocity  at  a  given  time  from  the  orbital  parameters. 
The  accelerations  are  computed  from  the  force  field  formulation.  The  inte¬ 
gration  is  performed  upon  the  derivatives  of  the  orbital  parameters. 

a.  Representation 

(1)  Initialization 


Prior  to  any  calculations,  the  epoch  elements  must  be  deter¬ 
mined  from  the  input.  Since  the  input  is  either  in  the  form  of  7  element 
cards  or  from  the  SEAI  tape,  the  system  subroutine  NXTELIi  is  used. 


performed: 


Given  the  epoch  parameter  set,  the  following  computations  are 


p  =  h  ■  h 
o  — o  — o 


W  = 


h 
— o 


—o 

•fh 

cos  i 

=  V 

0 

z 
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sin  i 

0 

=J~ 1  -  W, 

i 

=  arc  tan 

sin  i 


COS  1 


sinfl  =  — : — — 
o  sm  i 


cos  Q,  - 


o  sm  l 
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sm 


Q,  =  ate  tan  - — 

o  cos  Q 


2  2  ,  2 
e  -  a  +  a  „ 
o  xN  yN 

o  J  o 


a.  = 
o 


1  -  e 


U  =  L  -Q 
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if  W  >0 
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L  +  Q 
o  o 


if  <  0 
"'o 


•fir 


o  3/2 

a 
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a  =  a  „  cos  Q,  -  a  „  sin  cos  i 
x  xN  o  yN  o  o 

o  o  o 


a  =  a  „  sin  U,  +  a  „  cosQ  cos  i 
y  xN  o  yN  o  o 

-'o  o  J  o 


a  =  a  „  sin  i 
z  yN  o 

o  o 


q  =  a  (1  -  e  ) 
o  c  o 


If  the  input  mass  parameter  =  0,  mass  is  set  equal  to  1.  The  quantity  1/mass 
is  formed. 

(2)  Position  and  Velocity 

At  each  time  point,  whether  at  an  integration  step,  observation 
time  or  prediction  point,  a  similar  initialization  takes  place.  Then,  the 
program  goes  on  to  compute  the  geocentric  position  (r)  and  velocity  (r)  of 
the  satellite  .  0) 

Given  a,  h,  L,  compute  jt,  r_ 

p  =  h_  *  h_ 

2 

e  =  a  •  a 
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Kepler's  equation  is  solved  by  iteration  using  the  Newton-Raphson 
method  with  an  initial  guess  for  (E+UJ )  of  U  [i.e.,  (E+cn)1  =  U  ] 


[U  +  e  sin  E.  -  (E+co).  ] 

(E+C0)i+1  =  (E+CO). - e  cos'l" 1 - —  radians 

where 

e  cos  E.  =  a  lT  cos  (E+Oi).  +  a  ,,  sin  (E+a>). 

1  xN  1  yN  'r 

e  sin  E.  =  a  „  sin  (E+<u).  -  a  „  cos  (E+cn). 

1  xN  l  yN  v  'x 


The  iteration  is  concluded  when 
(E+03)i+i  -  (E-JO)  ).  <  10"8 

If,  after  50  iterations,  the  criterion  is  not  met,  the  run  is 
terminated.  A  comment  to  this  affect  is  written  on  the  output  tape. 

After  Kepler's  equation  is  solved,  the  calculations  continue: 

r  =  a  (1  -  e  cos  E) 


a  r  .  .  e  sin  E 

cos  u  =  -  [cos  (E  +cu)  -  a  +  a  - - 1  ' "n 

r  XN  yN  (1  +J 1  -  e  ; 


sin  u  =  ~  [sin  (E  +cn  )  -  j 

U  =  N  cos  u  +  H  sin  i 
V  =  -  N  sin  u  +  M  co. 
r  =  r  U 


e  sin  E 


u 


r  =  r  U  +  rv  V 
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b.  Drag 

To  compute  the  drag  perturbation,  it  is  first  necessary  to  calculate 
H,  the  altitude  above  an  oblate  spheroid: 

H  =  (r  -  1)  -  |  f2  (U  )4  +  (f  +  |  f2)  U  2 
Z  z  JL  z 


For  this  altitude  the  1962  U.S.  Standard  Atmosphere  Table  is  searched 
nsity  (p^)  and  molecular  weight  (Mg). 

Values  above  700  km  are  obtained  by  linear  extrapolation. 


The  acceleration  is  obtained  by  the  following  computation'1  ' : 
•  • 

V  =  x  +  y 6  where  9  =  0.058  834  47  radians/k  min 

x  e 

V  =  y  -  x@ 

y 

v2  =  * 


/  2  2 

2 

+  V 

x  y 

z 

.  P*  S 

^  4  ea 

s 

3  3 

V 

6.972  x  109 

vd 

,  1/4 


CD  **  MeTs 

O 


if  Cp  >  1418, 


otherwise, 


CD  =  CD  (1  +  Xl 1739130  e  ) 

0  2 
,  K  7T  d  , 


°c  =  CD  PrV  (- 
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The  acceleration  (rg)  due  to  the  zonal  harmonics  of  tae  Earth  s 
gravitational  potential  are  calculated  by 


S  Jn  (f>"p,„+l  <"*>] 

r  n=2 


V  y  t  3  *  Jn  (f>  p'n+l  (UZ>  ] 

r  n=2 


=  t  ■‘j  Jn  <f>“PVl  <".>) 

r  n=2 


-T  2  Jn(f)n^(V 

r  n=2 


»her8  P0  (02)  •  1,  Pj  «l2)  •  U2 


P„  (02;  *  ;  I  <2»  -  1)  %  pn-!  «„>  ‘  (Q’1)  'n-2  <V  1  ' 

\ 

p'  <V  ■  4“  "  (  “2  P„  <V  '  Pn-1  (V  ) 


The  coupling  factors  used  are 


d.  Radiation  Pressure 


The  acceleration  due  to  direct  solar  radiation  pressure.  r^, 

is  formulated  as  detailed  in  Reference  2.  Initially,  the  mean  longitude 
of  the  sun  at  the  beginning  of  the  current  year  L  is  found  by  the  system 
subroutine  TLC^); 
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If  cos  ip  >  0,  the  satellite  is  illuminated. 

I  £  cos  ip  <  0,  form  sin  (ip  +r))  =  sin  ip  cos  rj  +  sin  T)  cos  ip 


J {cos2  ip-  l)(-j  -  1)  +  (i)  cos  ip 


If  sin  Qp  +T))  >  0,  the  satellite  is  illuminated. 
If  the  satellite  is  illuminated,  calculate 
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e.  Derivatives  and  Integration 

The  total  perturbative  acceleration  can  now  be  used  to  determine  the 
perturbative  derivatives  of  the  parameters. 
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3.3  WEIGHTED  DIFFERENTIAL  CORRECTION 


It  is  necessary  to  simulate  observations  in  two  phases  of  the  SPWDC. 
The  most  obvious  phase  is  that  of  prediction.  Also,  in  the  differential 
correction,  the  estimate  of  the  observed  quantities  on  the  basis  of  the 
current  parameter  estimates  must  be  obtained.  This  estimate  is  often  referred 
to  as  the  "computed"  quantity  as  distinguished  from  the  "observed"  quantity. 
The  difference  between  these  quantities  is  the  "residual,"  e.g.,  for  the  ith 
range  observation 

A =  p (observed)  -  p.  (computed) 


Each  residual  is  related  to  corrections  in  the  parameters  by  an 
"equation  of  condition,"  e.g., 

Ap  =  2  AP. 

j=l  SPj  J 

where  p.  are  the  parameters  n^,  a  ,  etc. 

J  0  xw0 

These  equations  of  condition  are  weighted,  according  to  the  numbers 
obtained  as  outlined  in  Section  3.1,  and  processed  as  described  in  the  present 
section  by  relating  the  "computed"  quantities  to  the  geocentric  ephemeris 
obtained  as  described  in  Section  3.2.  This  relationship  involves  the  location 
of  the  center  of  the  Earth  with  respect  to  the  observing  station,  R: 

p  =  r  +  R 

The  partial  derivatives  with  respect  to  the  seven  parameters  are 
computed  according  to  analytical  formulas,  given  in  the  following  subsections. 
When  the  eccentricity  exceeds  a  given  maximum,  however,  the  coefficients  of 
„he  ballistic  parameter  are  computed  by  a  variant  ephemeris.  During  this 
calculation  the  quantities  p  ,  L  ,  and  p  are  stored  in  a  bloc  designated 
VDATA.  The  v  subscript  is  used  hereafter  \o  identify  these  results  of  the 
variant  calculation.  The  variant  ephemeris  differs  from  the  nominal  ephemeris 
by  using  a  ballistic  parameter  differing  from  the  nominal  value.  This  varia- 

/VI 

tion  in  the  ballistic  parameter  0 -  is  represented  in  the  following 
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r  1 

formulation  by  its  factor  0  (-) .  Therefore,  the  factor  C^A  appears  in  the 
coefficients.  The  correction  actually  should  be  thought  of  as  applying  to  not 
only  the  whole  ballistic  parameter  but  including  a  correction  to  the  mean 
density  as  well.  The  correlations  between  any  separate  determination  of  these 
factors  would  be  very  high  indeed  when  only  satellite  position  and  velocity 
observations  are  used. 

a .  Preliminary  Calculations 

This  part  of  the  program  calculates  "observations"  rf  the  orbit  spec¬ 
ified  by  the  input  set  of  elements  (or  subsequently  corrected  elements)  made 
from  the  sensor  station.  Given  the  station  coordinates  of  latitude,  6,  east 
longitude,  Xg,  and  height  above  sea  level,  H,  and  Greenwich  sidereal  time, 

G  ,  at  satellite  epoch,  t  ,  the  following  procedure  computes  the  observations 
to  •  o 

p  ,  p  and  the  direction  cosines  of  the  unit  vector,  L,  from  the  station  to  the 
c  c  — 

satellite  for  time,  t,  in  minutes  since  epoch. 

Compute  local  sidereal  time,  9,  at  time,  t: 

9  =  X  +  Q  +  (0.0043752691)  t  (Hod  2tt  ) 

h  to 

Compute  the  station  vector,  R: 


X  = 


r_JL_ 

L  cos  9 


] 


cos  9 


Y=  [ 


Y 

sin  9 


] 


sin  9 


Z  =  Z 

Compute  the  slant  range,  Pc: 
£c  =  r  +  R 

Pc=VPx  +  Py  +  Pz 
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p  =  x  +  X 

X 

p  =  y  +  Y 

y 

p  =  z  +  z 

z 

where  x,  y,  z  are  components  of  r_  obtained  from  Section  3.2  by  means  of 
input  elements  and  the  time,  t. 

Compute  unit  vector  from  the  station  to  the  satellite  in  the  equatorial 
coordinate  system: 


xc 


L2E 
Pc 

=  *L 


yc 


Lzc  P 


*1 

«-/ 
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Compute  range  rate,  p  , 

c 


X  =  -Y  9,  0  =  0.058,834,47 

Y  =  X  9 
Z  =  0 

4  =  I  +  R. 

4  =  L  Pc  =  Lx  (x  +  X)  +  Ly  (y  +  Y)  +  L^z, 
where  x,  y,  z  are  components  of  £,  obtained  from  Section  3.2. 
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b .  Range  Observations 


If  range  is  observed,  the  residual  is 
R1  Pc 

c.  Angular  Observations 

If  azimuth,  A,and  elevation,  h,  are  observed,  the  residuals  are, 
respectively, 

&2=  PCA*  (L-Lc) 

r3  =  Pc  D  •  <L  -  Lc) 

where 
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The  j>,  E,  Z  unii.  vector  system  and  the  horizon  oriented  L^,  A^,  unit  vector 
system  are  defined  by: 

S  =  sin  i  cos  0 
x 

S  =  sin  i  sin  9 

y 

S  =  -  cos  i 
z 
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If  right  ascension, &  ,  and  declination,  6  ,  are  observed,  the  resi- 
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duals  are: 
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L  =  cos  6  cos  Cl 
x 

L  /  L  =  cos  5  sin  a 

-  \  y 

L  =  sin  6 
z 

See  Figure  4  for  these  vector  relationships, 
d*  Range  Rate  Observations 

If  range  rate,  p  ,  is  observed, 

R6  =  pc  A p=  (  p  -  Pc)  Pc 
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FIG.  4  VECTOR  RELATIONSHIPS 
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e .  Differential  Correction  Solution 

The  residuals  in  the  observations  are  related  to  the  changes 
in  the  elements  through  a  first  order  approximation  of  the  form 


where  the  form  or  the  coefficients,  C^,  depend  on  the  observation  type  in 
residual,  time  of  observation,  and  the  observation  weights.  These 
coefficients,  which  have  been  developed  by  means  of  first-order  partials, 
are  functions  of  the  orbit  elements  and  computed  observations. 

The  coefficients  are  computed  by  first  establishing  the  R  and  U 
coefficients  at  time  t. 

2 

R  =  (a  /r)  e  sin  E 
u 

Rn  -  -2/3  r  +  (U  -  U0>  R, 

<a2/r>  [\,li  -  <E  +a>)] 

v (a2/r)  [> '  si° <E+a,>] 

U  -  (a2/r) 
u 


U 


n 


(U  -  U  )  U 
o  u 
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are  given 


the  range 
formula. 


xN  r 


(1  +  -)  sin  (E  +CO)  + 

a 


a  e  sin  E 
xN 


e*~  -  (1  +  l/l-e2)  e  cos  E 
(1  +  l/l-e2)2/'/l-e2 


yN 


1  W  1-e2 


UyN  ~  f  (1  +  cos  (E  +  a>)  + 


a  „  e  sin  E 
yN 


2  n 

e  -  (1  -!  Vl-e  )  e  cos  E 


xN 


1  +  V  1-e 2 


When  the  slant  range  is  observed,  the  residual  and  coefficients 
by  the  expressions: 


R  =  ~ 

°P 


C.  =fL.UR  +  1  •  V  u  "I 
An  1_— c  -  n  -c  —  n  J 

n 

CA„  =Tl.UR  +L.VU  1-4- 

AUq  J_— c  —  u  -c  —  u  J  °p 

=  [ic*  —  r  cos  ■*"  “  —  r  s^n  *  cos  U^J  "5" 


CA. 

Ai 


CL  •  W  r  sir.  ul 

-c  -  J  °p 


Ai 


Pv-Pc 

»  L6(;)J 


-0—  where  the  v  subscript  refers  to 


computed  from  the  variant  ephemeris.  If  computed  by  the  analytic 


’a| 


-  a  n  (t  -  t  )  C_  An  C,  K 
3  0  D  rr  An 
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where  K  is  a  unit  conversion  constant  and  is  set  e'U’al  to  2.  The 
quantity  00  m  the  coefficient  equations  is  the  standard  deviation  in 

the  range  measurement,  so  that  -g—  is  the  weight  of  that  measurement, 

p 

assigned  as  described  by  Section  3.1. 

R2 

When  azimuth  A  is  observed,  then  R  =  ~  and  the  coefficients 

A  1 

are  obtained  by  replacing  L  by  A.  In  this  case,  the  -  is  replaced  by 

'  ~  —  a 

1  P 

where  0  is  the  standard  deviation  in  the  azimuth  measurement. 

P  °,  A 

c  A 

The  seventh  coefficient,  in  the  variant  ephemeris  mode,  must  be  computed 
as 


.  fp  A*(L  -  L) 
1  rc  —  '-v  — c 

A  -  Pc  °A  <5  (!) 


With  the  azimuth,  elevation  h  is  always  observed,  and  in  this  case 
R 

R  =  —  and  the  corresponding  coefficients  are  obtained  by  replacing  L 

1  1 

with  D.  -q-  is  replaced  by  -p — g-,  where  0^  is  the  standard  deviation  in 

P  c  h 

the  elevation  measurement.  The  seventh  coefficient  in  the  variant  ephem¬ 
eris  mode  is 


P  °u 
rc  h 


p  D-(L  -  L  ) 
'c  —  — V  — c 


If  right  ascension  a  is  observed,  then  R  =  —  and  the  corres¬ 


ponding  coef f icie..“s  are  obtained  by  replacing  L  with  A.  -g —  is 
1  P 

replaced  by  p  •  —  ,  where  Oa  is  the  standard  deviation  in  the  right 

ascension,  ifte  seventh  coefficient  in  the  variant  ephemeris  mode  is 
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A± 

m 


p  O 
rc  a 


p  A  .  (L  -  L  ) 
c  -  — V  — c 


Along  with  right  ascension,  declination  6  is  always  observed,  and  in 
R5 

this  case  R  =  —  and  the  corresponding  coefficients  are  obtained  by 
replacing  L^,  with_D_.  -g—  is  replaced  by  p — ,  where  Og  is  the  standard 
deviation  in  the  declination.  The  seventh  coefficient  in  the  variant 
ephemeris  mode  is 


1 


A- 

m 


p  D  .  (L  -  L  ) 
*C  _  y  C/ 


6<s> 


If  range  rate,  p,  is  observed,  the  preliminary  coefficients 
are  computed  from: 


5/2  2 

R  ='/p  3  -  (e  cos  E  -  e  ) 
u  3 

r 


R  =  \  +  (U  -  U  )  R 

n  3  o  u 


5/2 


xN 


=«/>  ^-r-  [sin  (E  4co)  -  a  e  sin  E  -  a  ] 
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*yN  =  ^V"  [-  cos  f-E  +a>)  '  ayN  e  sin  E  +  axN  J 


Uu  "  — 

r 


-  e  e  sin  E 


U  =  ^  +  (U  -  U  )  U 
n  3  o  u 


5/2  / - r  r  /  2 

UXN  v 1  -  6  cos  (E  +CU)  -  axN  i1  +  rP 


U..M  *'frST‘  \/l  -  e2  Sin  (F.  +  to)  -  a  „  (l  + 


ap  / 


The  range  rate  coefficients  are  then  computed  as 

CA  =IL.U  ro  (R  -  vU  )  -  p  R  1+n.UR 
A  n  ^  -c  —  l  Hz  n  n  rc  n J  Cc  -  n 

n  ' 

+  [a  <u„  +  f  u„>  -K 

*  -1  V  2.  [  Pc  <**1.  •  *  V  ■  K  Rx»]+  K‘  £  “*N 

xN  [ 

+  L  ♦  V  |  p  (U  H  +  -  U  -  p  »J+P  'Vui^- 
— c  —  l *c  xN  r  xN  c  xN  J  — c  —  xN  I  up 

• 

CA  =  J  L  •  U  r  p  (R.-VU  -p  R  nl  +  P  •  UR  „ 

Aa  „  1  -c  —  l  rc  yN  yn  rc  yNJ  — c  —  yN 

yN  l 

+  L  •  V  rp  (U  „  +  -  U  ..)  -  p  UMl-tfi.VUu  ~ 

— c  —  L’c  yN  r  yN  c  yNJ  —  yN  >  Op 

CAI,  =J  L  *  U  [p  (R  -  v  U  )  -  p  R  l  +  p  -  UR 

AU  i  — c  —  lyc  u  u  <c  uJ  -c  —  u 

0  \ 


+  L  «  V  rp  (U  +  “  U  )  -  o  U  +n*Vul  — 4- 
-c  -  Inc  u  r  u  ^c  u-  Pc  1  ufdp 
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CAft  =  \  ~  Pc  -c-  r  *  C0S  1  +  V  -  cos  1  [  Pc  i  -  Pc  c  ' 

+  P ,»V  r  cos  i  +  L  «  W  sin  i  [p  (r  v  sin  u  -  r  cos  u 
— c  -  c  —  c 

+  P  rcosu]  -  n  >W  r  sin  i  cos  u[ 
c  -=c  -  I  CFp 

c  W  [Pc  (r  v  cos  u  +  r  sin  u)  -  r  sin  u  ] 

+  P  ■ 


■c  -  Op 

i  W  l  i 

s'  « 6  r<7> 


If  computed  by  the  analytic  formula, 


'si'1?*"*'-  '»>  CDA<>r  CA„KlST 


Let 


N 

X  C.  .  A.  =  R.  . 
j  =  l  J 


represent  all  such  equations  of  condition,  where  are  the  coefficients, 

R  are  the  accepted  observation  residuals,  A.  are  the  corrections  to  the 
ij  J 

orbital  elements,  and  N  is  the  number  of  elements  to  be  corrected  and  is  the 
number  of  accepted  observation  residuals.  The  resulting  matrix  equation  is 
solved  to  give  the  corrections,  A.  in  a  least  square  sense,  to  the  orbital 


elements  at  time. 

,  t 

.  These 

denote  corrected 

elements) : 

An 

n'  = 

II 

(1  +  “ 

0 

0 

n 

• 

U'  = 

0 

U 

0 

+  AU 

0 
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a 1 

=  a  . 

+  Aa 

xNo 

o 

* 

xN 

a1 

=  a  M 

+  Aa  „ 

yNo 

yNo 

yN 

tt'o 

=  ^o  + 

A  a 

i' 

0 

=  i  + 
0 

Ai 

L' 

=  u'  +n '  if  w 

0 

0 

0 

L* 

0 

=  U'  - 
0 

n;  if  w: 

i>  0 

i  <  0 


Following  the  above  calculation  of  the  corrected  elements, 
another  representation  of  the  observations  is  performed,  on  the  basis 
of  the  corrected  elements,  and  another  set  of  residuals  is  formed  by 
using  the  same  input  observations.  The  RMS  values  of  the  last  two  sets 
of  consecutive  residuals  are  compared  to  insure  convergence  of  the 
computational  process.  The  process  is  complete  when  the  residual  RMS 
converges  to  the  minimum  value  considered  as  acceptable. 
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3.4  PREDICTION  AND  RELIABILITY 


Once  a  corrected  set  of  elements  has  been  established,  the  Special 
Perturbations  formulation  is  used  to  compute  future  position  and  velocity. 

The  accuracy  of  the  prediction  may  be  estimated  by  means  of  an  analysis  of  the 
effect  of  the  element  uncertainties  by  their  propagation  through  the  represen¬ 
tation  equations. 

a .  Prediction  Options 


The  prediction  function  applies  the  equations  of  Section  3.2  to  repre¬ 
sent  the  position  and  velocity  at  some  specified  time.  Two  modes  of  prediction 
have  been  formulated:  (1)  prediction  by  specifying  the  initial  time, 
the  number  of  points  and  some  time  interval  and  (2)  prediction  for  the  time  of 
closest  approach  to  some  specified  sensor. 

Prediction  by  time  is  a  process  of  representing  the  position  and 
velocity  at  the  requested  times.  In  the  Special  Perturbations  formulation  this 
requires  an  integration  to  the  initial  time  point  followed  by  integration  to 
the  subsequent  time  points,  as  specified  by  the  time  interval,  until  the  final 
time  is  reached. 

Prediction  by  sensor  requires  a  more  complex  logic.  Here,  the 
numerical  integration  proceeds  until  the  range  rate  relative  to  the  station 
changes  sign  from  negative  to  positive.  At  this  point,  the  integration  is 
used  as  part  of  an  iterative  scheme  to  obtain  the  time  when  range  rate  is  zero. 
Once  established,  this  is  considered  to  be  the  time  of  closest  approach.  The 
routine  then  goes  on  to  obtain  prediction  points  for  specified  time  intervals 
and  duration  on  either  side  of  the  time  of  closest  approach.  This  process  may 
be  repeated  for  any  number  of  passes  over  the  sensor;  however,  it  is  restricted 
to  one  sensor . 

b .  Reliability  Estimates 

The  prediction  function  has  the  capability  of  estimating  the  reli¬ 
ability  of  its  prediction  points.  This  estimate  is  obtained  by  means  of  the 
variance-covariance  matrix  of  orbit  elements  established  in  the  final  iteration 
in  the  differential  correction  process  (Section  3.3  e) .  This  matrix  may  be  a 
six  by  six  or  a  seven  by  seven  depending  on  the  nature  of  the  differential 
correction.  In  either  case  the  variance-covariance  matrix  may  be  defined  as 

[A  *0  ,  and  the  quadratic  form 
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*7 

*- 


* 


m 


* 


«•* 


1 

1 


[D]  =  [  G]  [A*1]  [G]T 

leads  to  the  variance-covariance  matrix  of  predicted  position  [D], 

The  variance-covariance  matrix  of  predicted  velocity  [6]  is 
obtained  from  the  quadratic  form 


[D]  =  [H]  [A"1]  [H]T 


The  formulation  for  the  standard  deviation  in  predicted  position 
and  predicted  velocity  requires  the  formation  of  the  [G]  and  [H]  matrices.  The 
[g  ]  matrix  is  used  to  relate  uncertainties  in  orbital  elements  to  uncertainties 
in  predicted  position.  The  [h]  matrix  is  used  to  relate  uncertainties  in  orbital 
elements  to  uncertainties  in  predicted  velocity. 

The  [g]  matrix  is  composed  of  the  following  matrix  elements: 


S11  “ 


812  =  Rxn 

813  =  Ryn 


’14 


=  R 


317 
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g15  = 


816  =  0 


2/3  an(t-t  )  C.  A  p  R 
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g23  =  Uyn 
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u 

=  r  cos  i 
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o  D 

=  0 

=  0 

=  0 

=  0 
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n 


The  [h]  matrix  is  composed  of  the  following  matrix  elements 
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The  quantities  R  ,  R.  ,  R  ,  R  ,  U  ,  U  ,  U  ,  U  ,  R  ,  R  ,  R  , 
n  xn  yn  u  n  xn  yn  u  n  xn  yn 

R  ,  d  ,  U  ,  U  ,  U  are  defined  in  Section  3.3  e. 
u  u  xn  yn  u 

The  matrices  [ G  ]  and  [H  ]  are  functions  of  the  orbit  parameter 
and  the  prediction  time.  [  G  ]T  and  [H  f  are  the  matrix  transpose  of  [G  ] 
and{H],  respectively. 


Thus,  the  reliability  estimates  are  obtained  from  che  main  diag¬ 
onals  of  [D  ]and  [6].  These  estimates  are  the  standard  deviations  of  the 
in-track,  cross-track,  and  out-of-plane  components  of  position  and  velocity 
at  the  prediction  time. 
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SECTION  4 


OPERATING  INSTRUCTIONS 


The  Executive  Program  of  the  SPS  B-2  System  processes  the  Element, 
Sensor,  Observation  and  Parameter  cards  leaving  the  data  ir.  specified  B-2 
core  buffer  locations  to  be  accessed  by  OBSWGT  and  SPWDC. 

The  Observation  Weighting  Program  will  reformat  the  observations, 
ordering  them  in  time  and  applying  weights  (optional)  for  a  maximum  of  400 
observations.  It  then  calls  SPWDC  for  execution. 

The  SPWDC  program  may  go  through  a  differential  correction  process, 
prediction  by  station,  or  prediction  by  time  as  specified  on  the  Program 
Parameter  Cards. 


4.1  TAPE  SETUP 


The  logical  tape  assignments  are: 


Logical  Unit 


Tape  Description 


1 

2 

4 

6 

7 

11 

12 


SPS  B-2  Binary  Master 

Input 

SEAI 

SRADU  (optional  observation  input) 
Weight  tape  (optional) 

Output 

r,  i  binary  output  tape  (optional) 


4.2  DECK  SETUP 


Input  to  the  SPWDC  Program  consists  of  SPS  control  cards,  Parameter 
cards,  Element  cards,  Sensor  cards  and  Observation  cards  (Figure  5). 


FIGURE  5.  SPWDC  INPUT  DECK 


4.3  INPUT 

a.  Program  Notes 


(1)  The  input  options  for  the  SPSJOB  card  are  to  be  found  in 
Appendix  I, 

(2)  The  differential  correction  is  controlled  by  the  following 
cards: 

(a)  The  OBSWGT  Parameter  card  calls  the  SPWDC  Program. 

(b)  The  Vehicle  Characteristic  card  contains  the  mass, 
diameter,  and  reflectivity  of  the  satellite  for 
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the  drag  and  radiation  pressure  computations. 

This  card  may  be  omitted  if  the  following  values 
are  acceptable: 

mass  =  10.0  kg 

diameter  =  1.0  meter 

reflectivity  =  1.0 

(c)  The  Differential  Correction  Control  card  contains 
the  pertinent  information  for  the  integration  and 
differential  correction. 

(d)  The  Absolute  Error  Control  card  may  be  omitted  if 

aL  ■  l<f7 

V  \  ■  10’8 

are  acceptable  values. 

(e)  The  Relative  Error  Control  card  may  be  omitted  if 


are  acceptable  values. 

(f)  The  Program  Execution  card  specifies  that  a  differen¬ 
tial  correction  is  to  be  executed  and  starts  opera¬ 
tion  of  the  program. 

(3)  Prediction  by  time  is  controlled  by  the  following  cards: 

(a)  OBSWGT  Control  card 

(b)  Absolute  Error  Control  card 

(c)  Relative  Error  Control  card 

(d)  Time  Prediction  card  #1.  (This  card  contains 
integration  and  output  information.) 
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(e)  Time  Prediction  card  #2.  (The  time  output  interval 
is  indicated  on  this  card.) 

(f)  Program  Execution  card.  (The  Prediction-by-time 
path  is  set,  and  execution  begins.) 

(4)  Prediction-by-station  is  cc .'.trolled  by  these  cards: 

(a)  OBSWGT  Parameter  card 

(b)  Absolute  Error  Control  card 

(c)  Relative  Error  Control  card 

(d)  Station  Pass  Prediction  card  #1.  (The  integration 
and  output  information  is  specified.) 

(e)  Station  Pass  Prediction  card  #2.  (Sensor  data  may 
be  input  from  this  card  or  from  a  B-2  SEAI  tape.) 

(f)  Program  Execution  card.  (The  Prediction-by-sensor 
path  is  set,  and  execution  begins.) 

(5)  Each  of  the  three,  differential  correction,  prediction- 
by- station,  or  prediction-by-time,  may  be  run  separately, 
or  the  combination  of  differential  correction  and  pre¬ 
diction-by-station  or  time  may  be  executed.  For  actual 
test  cases  and  variations  of  logical  oaths,  see  Section  6. 

b .  Input  Deck 

A  typical  input  deck  is  arranged  as  follows: 

70  SCHTP  (U82) 

JOB  Card 
REM  Card 

SPSJOB 


Parameter  Cards 
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a 


\ 


M  f 


ir 

? 


i 


Element  Cards  (SPADATS  format) 

(3) 

Sensor  Cards  (SPADATS  format) v  1 

(2) 

Observation  Cards  (SPADATS  format) v  1 


ENDOFJOB 

(U82) 

ENDSCHED 

(U82) 

Since  the  Element,  Sensor  and  Observation  cards  are  standard 
SPADATS  System  cards,  only  the  Program  Parameter  cards  and  the  Weighting 
Tape  Input  will  be  described  in  this  Section. 

(1)  Parameter  Card  Formats 


The  first  parameter  card  must  be  the  OBSWGT  card  which 
contains  the  identity  of  the  program  called  by  OBSWGT  after  processing 
the  observations.  The  rest  of  the  parameter  cards  may  be  in  any  order, 
with  the  exception  of  the  Program  Execution  card  which  must  be  the 
last  card.  Note  that: 

(a)  If  six  or  seven  elements  are  being  corrected, 
the  n  only  correction  and  the  A<3  check  may  be 
chosen  as  options.  (Differential  Correction 
Control  card) 

(b)  The  prediction  reliance  option  for  prediction- 
by-time  or  station  is  only  possible  if  six  or 
seven  elements  are  to  be  differentially  corrected, 
the  observations  are  weighted,  and  the  drag  per¬ 
turbation  is  computed  (for  correction  of  seven 
elements) .  If  these  conditions  are  not  met,  the 
flag  is  turned  off  and  the  run  proceeds  without 
prediction  reliance.  (Time  and  Station  Prediction 
cards) 

Appendix  I  contains  the  Parameter  card  formats  and  related 
code  information. 


1 
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1 


(2)  Weighting  Cards 


Additional  weighting  input  data  may  be  contained  on  a 
weighting  tape  which  has  the  following  card  deck  setup: 


FIGURE  6 .  CARD  DECK  FORMAT  FOR  THE  WEIGHTING  TAPE 

Data  will  consist  of  the  standard  deviations  (1)  f6r  the 
observations  specified  as  constants  for  each  sensor,  (2)  for  each 
observation,  or  (3)  specified  by  control  information  for  determining 
the  method  to  be  used  for  computing  the  standard  deviation. 

There  may  be  weighting  of  all  observations  jor  no  weighting, 
but  in  a  particular  job  there  cannot  be  both  weighted  and  unweighted 
observations.  There  must  be  sigmas  for  all  quantities  observed  or 
the  observation  will  be  rejected  for  lack  of  weighting  information. 
Weighting  information  is  obtained  from  Weight  tape  or  from  the  file 
data  in  the  Observation  Weighting  Program.  (Appendix  I  contains  the 
card  formats  for  the  Weighting  tape.) 

A  weighting  tape  must  contain  weighting  information  for 
(1)  sensors  only,  (2)  observations  only,  (3)  functions  only, or  (4) 
sensors  and  functions. 
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The  output  for  SPWDC  consists  of  element  cards,  printed  data,  and 
optionally  a  binary  tape  and/or  punched  cards.  The  SPADATS  seven  element  cards 
plus  the  vehicle  characteristic  card  are  punched  out  (data  select  2,  code  mode) 
after  the  successful  completion  of  the  differential  correction  of  the  orbit,  A 
binary  tape  containing  time  in  minutes  since  epoch,  position  (r)  and  velocity 
(r)  is  an  optional  output  for  prediction.  See  Figure  7, 


18  sets  of  t,  r,  r 
followed  by  two  words 
of  D/0 


1-18  sets  of  t,  r,  r 
followed  by  two  words 
of  A/ZZZZZZZZ 


FIGURE  7.  CONTENTS  OF  BINARY  TAPE 

The  printed  data  (data  select  1)  consist  of  heading  lines,  differential 
correction  output,  prediction  output  and  error  comments.  An  additional  punched 
card  option  will  output  special  format  cards  containing  time  and  position  data. 
The  punched  cards  contain  £  and  integral  days  since  beginning  of  the  year  (first 
card)  and  fractional  time  (second  card)  in  the  format 

ID  +  XXXXXXXXX  +  YY 
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where  ID  is  an  identifier: 


f 

ID 

«-■ 

08 

X 

09 

Y 

? 

10 

Z 

jr 

11 

t 

(integral  days  since  beginning  of  year) 

12 

t 

(fractional  day) 

T 

I 

a.  Initial  Output 

*r> 

The  first  line 

of 

each  page  will  be: 

■1 

4- 

SPECIAL  PERTURBATIONS  WITH  WEIGHTED  DIFFERENTIAL 

CORRECTION  PRG  (Date)  PAGE_ 


This  is  followed  by  the  card  image  of  the  parameter  cards.  Any  error 
on  an  input  parameter  card  will  cause  the  following  comment  to  be  printed: 

ERROR  IN  INPUT  CARD  NO  JOB  TERMINATED 


M 


«r 


1 

ttjbv 


If  the  proper  sensor  data  was  not  found  on  the  SEAI  tape  or  on  a  sensor 
card,  the  remark 

NO  COORDINATES  FOR  STATION  NO _ ,  OBS.  SKIPPED 

is  printed.  When  no  standard  deviations  are  found  for  an  observation, 

NO  WEIGHTING  INFORMATION  FOR  OBS.  FROM  STATION 
NO _ ,  OBS.  SKIPPED 

is  printed. 

b .  Differential  Correction  Output 


1 

«** 


1 

1 


The  differential  correction  output  consists  of  the  residuals  (optionally) 
for  all  of  the  observations,  followed  by  the  root-mean-squares,  correction  to  the 
elements  and  corrected  elements.  For  a  weighted  run  the  normalized  root-mean- 
squares  are  also  output.  This  output  is  repeated  until  the  differential  correction 
converges  or  diverges . 
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For  convergence  the  finai  set  of  residuals  is  followed  by  the  comment 

DC  CONVERGED  -  THE  NEXT  CORRECTIONS  WOULD  BE  .  .  . 

and  the  final  set  of  corrections.  This  is  followed  by 

NO.  OF  RESIDUALS  USED  =  _ NO.  OF  RESIDUALS 

REJECTED  =  _ 

printed  with  the  values.  Then  follows 

OLD  ELEMENTS  WITH  RESPECT  TO  OLD  EPOCH  .  .  . 

OLD  ELEMENTS  WITH  RESPECT  TO  INTERMEDIATE  EPOCH  .  .  . 

NEW  ELEMENTS  WITH  RESPECT  TO  INTERMEDIATE  EPOCH  .  .  . 

NEW  ELEMENTS  WITH  RESPECT  TO  NEW  EPOCH  .  .  . 

with  appropriate  values.  The  seven  element  cards  and  the  vehicle  characteristic 
card  are  punched  out  when  the  differential  correction  has  converged.  Following 
this  output,  the  comment 

ELEMENT  SET  UPDATED 

is  printed.  For  divergence,  the  error  comment 

DIFFERENTIAL  CORRECTION  DID  NOT  CONVERGE,  JOB  TERMINATED 

is  followed  by  the  proposed  corrections  to  the  elements  and  the  corrected  elements 
with  respect  to  the  old  epoch. 

c.  Prediction  Output 

For  prediction  by  print  triples  there  are  four  output  options.  The 

time  in  days  since  beginning  of  the  year  is  followed  by  any  or  all  of  the 

following: 


(1) 

position  and  velocity  (r,  _r) 

» 

(2) 

osculating  orbital  elements 

(s,  6,  ij  £2  >  CD  j 

U) 

(3) 

subsatellite  track  (0,  X 

h) 
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{Uuwgk 


output 


For  prediction  by  station  pass  there  is,  in  addition,  the  option  to 


T 


•tfo 


* 


(4)  acquisition  coordinates  (  p, p ,  A,  h) 

d.  Prediction  Reliability  Output 

The  optional  prediction  reliability  output  will  print  the  standard 
deviation  in  the  predicted  position  and  velocity  components  at  each  prediction 
point.  The  format  provides  two  lines  of  data  per  time  point.  On  the  first  line 
the  standard  deviations  are  referred  to  the  cross-track,  in-track  and  out-of¬ 
plane  coordinates,  and  on  the  second  line  they  are  referred  to  the  inertial  x, 
y,  z  coordinate  system. 

e .  Error  Comments 


The  following  are  error  comments  which  may  appear  on  the  printed  output: 

(1)  An  error  in  the  Runge-Kutta  integration  routine 
will  give  one  of  the  following  comments: 

ERROR  IN  INTEGRATING  VARIANT  EPHEMERIS, 

JOB  TERMINATED.  (Correcting  Mass) 

ERROR  IN  INTEGRATING  EPHEMERIS  FOR  DC, 

JOB  TERMINATED.  (No  Mass  Correction) 

ERROR  IN  RUNGE-KUTTA  IN  PREDICTION  BY 
TIME  RUN  TERMINATED.  (Predict. on  by  Time) 

A (I)  +  R(I) Y(I)=Q  ERROR  IN  RUNGE  KUTTA 
INTEGRATION  (Prediction  by  Station) 

DELTA  T  =  0  AT  ANY  PT  IN  THE  INTEGRATION 
PROCESS  (Prediction  by  Station) 

(2)  In  computing  the  satellite  altitude  for  the  drag  per¬ 
turbation  for  too- low  a  satellite,  this  comment  is 
printed: 

SATELLITE  DROPPED  BELOW  50  KMS .  TIME 
SINCE  EPOCH  = _ MIN. 
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(3)  When  all  the  observations  are  rejected  before  the  differ¬ 
ential  correction, 

NO  GOOD  OBSERVATIONS  END  OF  JOB 

is  printed. 

(4)  If  Kepler's  equation  could  not  be  solved  for  (E  +cd) 
before  entering  the  integration  routine,  the  comment 

DID  NOT  CONVERGE  IN  KEPLERS  EQN  IK  50 
ITERATIONS,  RUN  TERMINATED 

appears  on  the  output. 

(5)  If  the  number  of  good  observations  was  less  than  the 
number  of  elements  to  be  corrected, 

DC  FAILED.  NOT  ENOUGH  GOOD  OBSERVED  QUANTITIES 
TO  SOLVE  THE  LEAST  SQUARES  MATRIX,  JOB  TERMINATED 

is  printed. 

(6)  If  the  RMS  gets  larger  twice  in  a  row,  if  there  is  a 
change  greater  than  0.5%  in  the  RMS,  or  if  the  maximum 
number  of  passes  through  the  differential  correction 
have  been  made,  the  following  is  printed: 

DIFFERENTIAL  CORRECTION  DID  NOT  CONVERGE, 

JOB  TERMINATED. 

(7)  When  the  input  Aq  is  greater  than  the  computed  Aq, 
this  comment  appears: 

DELTA  Q  TOO  LARGE,  REPEAT  CORRECTION  WITHOUT 
CORRECTION  AXN  AND  AYN. 

(8)  When  the  prediction  reliability  check  is  attempted  in 
a  correction  involving  less  than  six  elements,  this 
comment  appears: 

INSUFFICIENT  EQUATIONS  TO  OUTPUT  SIGMAS 
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SECTION  5 


FLOW  DIAGRAMS 


The  following  flow  diagrams  display  the  computational  procedures 
and  logical  flow  used  in  the  program.  Standardized  symbols  are  maintained 
throughout.  Ovals  are  used  to  indicate  subroutines.  Rectangles  are  used 
for  computational  processes.  Diamonds  represent  logical  decisions  or 
branching  tests. 


5.1  MAIN  FLOW  DIAGRAMS 

The  main  flow  of  the  program  is  shown  in  macro  structure  on  the 
following  two  pages.  These  diagrams  display  the  differential  correction 
function  and  the  prediction  function  of  SPWDC. 


5.2  DETAILED  FLOW  DIAGRAMS 

The  detailed  flow  diagrams  are  displayed  in  Appendix  II  and  show 
all  program  loops  within  SPWDC. 
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SECTION  6 


COMPILED  TEST  CASES 


The  three  test  cases  for  Satellite  116  described  in  this  section  will 
test  the  logical  paths  of  the  program.  Note  that  the  Observation  and  Element 
cards  are  the  same  for  all  cases  and  the  sensors  are  obtained  from  the  SEAI 
tape.  The  vehicle  characteristics  for  61  OMICRON  1  are  diameter  1,1283  m,  and 
mass  79.378  kg. 


6.1  TEST  CASE  1 

See  Figure  9,  Weighted  Differential  Correction. 

The  first  test  case  is  a  weighted  differential  correction  with  11 
observations  covering  a  40-hour  period.  The  element  set  and  observations  are 
card  inputs  to  the  program;  the  sensors  are  obtained  from  the  SEAI  tape  (input 
option  5).  The  weighting  data  are  from  the  Weight  Tape  (Figure  15).  Bulge, 
drag,  and  radiation  pressure  perturbations  are  included  in  the  variable  inte¬ 
gration  with  a  starting  step  size  of  -1.0  minute.  Six  elements  are  to  be 
corrected  in  a  maximum  of  9  passes  through  the  differential  correction  with  an 
n  only  correction  on  the  first  pass. 

The  absolute  maximums  for  rejection  are  1000  km  and  0.5  km/ sec;  the 
rms  multiplier  is  1.5.  All  residuals  are  to  be  output,  and  the  new  epoch  is 
the  10030  revolution. 

See  Figure  10  for  Test  Case  1  output. 
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6.2  TEST  CASE  2 


See  Figure  11,  Weighted  Differential  Correction  and  Prediction  by 

Time. 

The  ephemeris  calculation  incorporates  a  variable  negative  inte¬ 
gration  with  an  initial  step  size  of  one  minute  taking  into  account  bulge 
and  drag  perturbations  and  taking  the  new  epoch  to  be  200.7  days  since  the 
beginning  of  the  year.  Six  elements  are  to  be  corrected  with  an  n  only 
correction  on  the  first  pass  and  a  Aq  check  of  500  km.  The  time  prediction 
integration  is  in  the  variable  mode  with  an  initial  step  size  of  1.0  minute 
taking  into  account  bulge  and  drag  perturbations.  Ten  points  are  output  at 
a  one-minute  interval  starting  at  200./  days  (July  19,  1963).  All  printed 
options  including  prediction  reliance  values  are  output  as  well  as  the  binary 
tape  and  punched  cards. 

See  Figure  12  for  Test  Case  2  output. 


6.3  TEST  CASE  3 


See  Figure  13,  Unweighted  Differential  Correction  and  Prediction  by 

Station. 

This  is  an  unweighted  differential  correction  run  updating  the  epoch 
to  the  time  of  the  last  observation.  The  integration  is  in  the  variable  mode 
with  an  initial  step  size  of  -1.0  minute;  bulge  and  drag  perturbations  are 
included.  Six  elements  are  differentially  corrected  with  an  n  only  correction 
on  the  first  pass.  The  initial  rejection  maximums  are  1000  km  and  0.5  km/sec 
with  an  rms  multiplier  of  1.5.  There  is  a  maximum  of  nine  passes  through  the 
differential  correction  with  the  residuals  output  each  time. 


In  the  station  prediction  the  integration  mode  is  variable  with  a 
two-minute  step  size;  the  Earth's  bulge  and  drag  are  included.  Printed  output 
consists  of  t;  r,  r;  a,  e,  i,ft,co,  U;  0,  X  ,  h;  and  p,  p,  A,  h  for  sensor  216 
(coordinates  are  on  the  SEAI  tape).  The  output  points  will  be  5  minutes  on 
either  side  of  the  point  of  closest  approach  at  a  1.0  minute  interval  provided 
h  >1.0°. 


See  Figure  14  for  Test  Case  3  output. 
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FIGURE  10.  OUTPUT,  TEST  CASE  1 


FIGURE  11.  WEIGHTED  DIFFERENTIAL  CORRECTION  AND  PREDICTION  BY  TIME,  TEST  CASE  2 
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FIGURE  13.  UNWEIGHTED  DIFFERENTIAL  CORRECTION  AND  PREDICTION  BY  STATION,  TEST  CASE  3 
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FIGURE  14.  OUTPUT,  TEST  CASE  3 


WEIGHT  TAPE 


APPENDIX  I 


follows: 


PROGRAM  PARAMETER  CARDS 


Code  information  pertinent  to  the  parameter  card  formats  is  as 


Code 

1 

2 

3 


4 


Contents 


Alphanumeric  information 
Integers 

Floating  point  numbers 

XX. . .X.XX. . .X  (decimal  anywhere 
in  the  field) 

+.XXX. . .X+YY 

XXXXX . X  (no  decimal) 

Flag 

A  or  0:  not  to  be  computed 
1:  to  be  computed 
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3333333313$ 
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Hill 

HHIHI 
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Hill  Hill  Hill 

11111111111 
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17777777 
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lilllllllll 
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IIIHIII 
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IIIHIII 
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Column 

1-8 

9-16 

17 

18 


Contents 
SPSJOB 
OBSWGTE  A 

INPUT  OPTIONS  (See  following  page) 
0  OUTPUT  OPTION 


FIGURE  16.  SPSJOB  CARD 
-6-3- 


SCHEDULE  TAPE  INPUT  OPTIONS 


Sensor 

Option  Cards 


0 

1 

2 

3 

4 

5 

6 

7 

8 
9 


N 

N 

Y 

Y 
N 
N 

Y 

Y 
N 
N 


(Reference  SPSJ0B  Cards) 


Element  Observation 

Cards  Cards 


N 

Y 
N 

Y 
N 

Y 
N 

Y 

Y 
N 


N* 

N* 

N* 

N* 

Y 

Y 

Y 

Y 
N 
N 


Parameter 

Cards 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 

Y 


Satellite 

Cards 

Y 
N 

Y 
N 

Y 
N 

Y 
N 
N 

Y 


N  «=  No 
Y  =  Yes 


Options  0-7  refer  to  differential  correction  and  (optionally)  prediction 
8-9  refer  to  prediction  only 

♦Observations  are  from  the  SRADU  tape. 


NOTE:  Whenever  N  appears,  data  are  obtained  from  data  files  (except  for 
satellite  number  cards).  In  the  case  of  sensors  and  elements,  the 
data  source  is  the  SPS  B-2  SEAI  tape.  When  elements  are  input  from 
the  SEAI  tape,  the  Schedule  Tape  input  must  contain  satellite  number 
cards.  If  observation  cards  are  not  input  from  the  schedule  tape, 
they  are  obtained  from  the  SRADU  tape. 
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Field 

Column 

Contents 

1 

1  -  8 

SPWDC  AAA 

2 

9-11 

Satellite  Number  (optional) 

3 

12 

0  or  A  -  Weighting  Tape  Input 

1  -  Use  Sensor  Weighting  Data  in  OBSWGT 

Weight  File 

4 

13  -  78 

Not  used 

5 

79 

0  or  A  -  Use  Observation  Weights  in  SPWDC 

1  -  -o  Not  Use  Observation  Weights 

6 

80 

P 

FIGURE  I.  .  OBSWGT  PARAMETER  CARD 
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Code 
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1  - 
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52 

Diameter  of  Satellite  (Meters) 
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3 

53  - 
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Mass  of  Satellite  (Kg) 
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7  (Reflectivity) 
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FIGURE  18.  VEHICLE  CHARACTERISTIC  CARD 
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At  Integration  Mode  (A-  variable,  1-fixed) 
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Bulge  Perturbation 
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Drag  Perturbation 
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Radiation  Pressure  Perturbation 
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13 
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15 
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n  (RMS  Multiplier)  First 

Residual  Output  (A  -Never) (1  -  Set) 

lidS  L 

(2-Every  Set) 

3 

16 

72 

17 

73  -  77 

Not  Used 

18 

78  -  79 

A  3 

19 

80 

P 

*Applies  only  if  6-7 

Elements  are  to  be  corrected. 

FIGURE  IS.  DIFFERENTIAL  CORRECTION  CONTROL  CARD 
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FIGURE  2; .  ABSOLUTE  ERROR  CONTROL  CARD 
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FIGURE  21.  RELATIVE  ERROR  CONTROL  CARD 
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Field 

Column 

Contents 

Code 
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1 

At  Integration  Mode  (A-  variable) 

(1  -  fixed) 

2 

2 

2  - 

11 

At  (Min) (may  be  negative  if  desired) 

3 

3 

12 

Bulge  Perturbation 

4 

4 

13 

Drag  Perturbation 

4 

5 

14 

Radiation  Pressure  Perturbation 

4 

6 

15  - 

16 

Print  Flag  1-t,  it,  r_  . 

<  2-t,a,e,i L 
'  4-t  ,$ ,  Xg,  h 

Add  flags  for 
combination 
of  data 

7 

17 

Binary  Tape  Output 

4 

8 

18 

Prediction  Reliability* 

4 

9 

19 

Punched  Cards  (x_, t) 

4 

10 

20  - 

77 

Not  Used 

11 

78  - 

79 

A6 

12 

80 

P 

*Must  be  a  weighted  6  element  or  7  element  with  drag  differential 
correction  ruu. 

FIGURE  7.2.  TIME  PREDICTION  CARD  #1 
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Column 
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Code 
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1  - 

12 

Time  in  days  since 
beginning  of  year 

3 

2 

13  - 

20 

A  t  (Min) 

3 

3 

21  - 

24 

Number  of  output  points 

3 

4 

25  - 

77 

Not  Used 

5 

78  - 

79 

A  7 

6 

80 

P 

FIGURE  2"\  TIME  PREDICTION  CARD  #2 
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Column 

Contents 

Code 
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1 

At  Integration  Mode  ^  fixed)^ 

2 

2 

2  - 

11 

At  (Min) 

3 

3 
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Bulge  Perturbation 

4 
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13 

Drag  Perturbation 
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14 

Radiation  Pressure  Perturbation 
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Add  flags  for 
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of  data 

8-t,  p,p  ,  A,  h 

7 

17 

Binary  Tape  Output 

4 

8 

18 

Prediction  Reliance* 

4 
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Card  2  needed) 
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1 

1.C 

23  - 

32 

K  -  No.  of  Passes 

2 

11 

33  - 

42 

AT  (Min) (Either  Side  of  Closest  Approach)  3 

12 

43  - 

52 

At*  (Min)  (Output  per  At*) 

3 

13 

53  - 

62 

Minimum  h  (Deg) 

3 

14 

63  - 

77 

Not  Used 

15 

78  - 

79 

A  8 

16 

80 

P 

*Must  be  a  weighted  6  element  or  7  element  with  drag  diffeiential 
correction  run. 


FIGURE  24  .  STATION  PASS  PREDICTION  CARD  #1 
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FIGURE  25 .  STATION  PASS  PREDICTION  CARD  #2 
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APPENDIX  III 


SYMBOL 

PROGRAM  GLOSSARY* 

DESCRIPTION 

A 

a 

ACCCNT 

Total  residual  count 

AE 

a 

ALFLG 

e 

A,  a  residual  flag 

APRINT 

Output  location  for  a0 
d2 

Satellite  area  (n“) 

AREA 

ASTK 

* 

AXGR 

V 

a 

AXO 

X 

a 

x 

AYGR 

a'° 

AYO 

y 

ay 

AZGR 

a'0 

AZIM 

z 

Azimuth 

AZO 

a 

BFLAG 

z 

0 

Bulge  flag 

BGOB 

Reformatted  0BL0C  (16  word  format) 

BTFLAG 

Binary  tape  flag 

CAPD 

D  (DERIV) 

CC 

C  (DRAG) 

GGC40 

Switch  for  CCCC  first  time 

*Includes  (1)  SPWDC  program  symbols  and  (2)  all  symbols  not  appearing 
in  che  SPS  B-2  Assign  Deck.  Underscored,  capitalized  notations  refer 
to  subroutines  in  the  program. 
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APPENDIX  C  (Continued) 


CC.C50 

CDO 

CD 

CFLAGI 

CFLAG2 

C FLAG 3 

CFLAG4 

CFLAG5 

CNFLAG 

CNTRLX 

CONTEST 

COSEPS 

COS  PS  I 

COUNTLR 

COUNTL 

COUNTR 

DC  APT 

DCFLAG1 

DCONTI 


DCONT2 

DC 

DDGR 
DELTAO 
DELTATO 
DEN  7 


Switch  for  CCCC  n  only  correction  on  first  pass 

CD  =  0.92  '.DRAG) 

cd°(drag) 

First  time  through  flag  (CNTRL) 

No  more  observations  flag  (CNTRL) 

Variant,  normal  ephemeris  flag  ( CNTRL) 

New  epoch  by  revolution  or  time  flag  ( CNTRL) 

New  epoch  updated  flag  (CNTRL) 

Correct  n  flag 

Switch  for  D.  C.  Control  Region  (CNTRL) 
Convergence  rest  (5%)  (CCCC) 

cos  e  =  0 .9174469 
cos  ip 

Same  as  COUNTL 

Number  of  elements  being  corrected 
Same  as  COUNTL 


A  T  (Station  Prediction) 


Differential  Correction  flag 
„  3 


co 
4  ea 


3,  .  3  o  4 

cm  kemm  K 


gm  radii" 


=  2.419012831x10 


21 


-0.6378150x10 


10 


D 


3  .  ‘ 
cm  kg 


gm  radii 


D  ( DERIV) 
Aq 
A  t 
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DENOM 

1 

DFLAG 

DGR 

d 

DIAM 

fr 

DIVFL 

d 

DLFLG 

if 

DLT1M 

7 

DLVX 

DLVY 

DLVZ 

*- 

DQFLAG 

V  - 

DQN 

• 

DRSDL 

DTASK 

- 

DTDCRS 

DT 

— 

DTSROPT 

* 

E2VGR 

* 

J 

ECSIG 

*» 

ELEV 

1 

«*» 

ELFLAG1 

ELFLAG 

1 

EOCHK 

1 

1 

1 


-  3 


I  o 

1  +  Vl  -  e  (SUBXYZ) 

Drag  flag 
D'  (DERIV) 

Effective  diameter  of  the  satellite 
Divergence  flag 
6  ,  h  residuals  flag 
A  (^)  (DC0RR2) 

AVx 

AVy 

AVZ 

Delta  q  check  flag 
Maximum  A  q  for  A  q  check 
P  D*AL  (kms) 

•k 

A  t  (Station  Prediction) 

A  t  decrease  flag  (Station  Prediction) 

At  (prediction  by  time) 

A  t  option  flag  (Station  Prediction) 

2  v' 

-e  V 

CO 

e 

h  (RHODAH) 

Error  flag  for  prediction  reliance 
Prediction  Reliance  Flag 

Eccentricity  check  -  variant  ephemeris  vs.  analytical 
computation  for  mass 
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APPENDIX  C  (Continued) 


EPHEM 

EPRINT 

EPS1 

FOBHR 

FOBMIN 

FOBSEC 

FRDAY 

FRF LAG 

FRF RAC 

FRSTFLG 

FRTIME 

FSTCTLS 

FSTDTSR 

FSTPTW 

FSTSOUT 

GAMMA 

H2MH1 

H2MH 

HCONI 

HMH1 


Buffer  for  saving  ephemeris  L,  a,  h 
Output  location  for  eQ 
e  (RAD)  for  Kepler's  Eqn  =  10*8 

Output  location  for  observation  time  -  hr 

Output  location  for  observation  time  -  min 

Output  location  for  observation  time  -  sec 

Integral  days  since  beginning  of  year 

Punched  card  output  flag  (SUBOUT) 

Fractional  day 

First  time  flag  (Station  Prediction) 

Days  since  beginning  of  year  (SUBOUT) 

First  time  flag  (Station  Prediction) 

First  A  t*  Pass  flag  (Station  Prediction) 

First  printed  output  flag 

First  rrinted  output  flag  (Station  Prediction) 

Reflectivity  of  satellite 

H2  -  Hj  (TBINT) 

H2  -  H  (TBINT) 

3/2  f2  =  0.0168571736  x  10-3  (CALH) 

H  -  Hx  (TBINT) 
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HMIN 

HSUBQP 

INELT2 

INELT 

INP20 

IPRINT 

JBUF 

KEORTM 

KNTRL 

LINECT 

LPRINT 

LSQBUF 

L3TT 

LTBUF 

MASSO 
MASS PR 
MASS 
MAX  2 


h  .  (Station  Prediction) 
mi«i 

Output  location  for  HqQ  (miles) 

Beginning  storage  location  for  original  elements 
Beginning  storage  locations  for  updated  elements 
Switch  for  Parameter  Card  Input  (INPUT) 

Output  location  for  i  (deg.) 

Kozai  Values  for  J ^  through 

K e/Jjl  =  0.07436574 

Which  elements  to  connect  (m,  i ,  ft ,  U0’ayN’  axN’  n') 
Line  Count 

Output  location  for  (deg.) 

Least  Squares  Buffer 
Last  t  (PR00BS) 

Buffer  for  Controls  for  prediction 

Original  mass  of  satellite  (Kgs) 

Output  location  for  mass 
Program  value  of  mass 
p  rejection  criteria 
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MAX 

ME 

MFLAGl 

MSPRK 

MU 

NLINE 

NOPASS 

NPRINT 

N 

NH 

NUX 

NUY 

NUZ 

OBDAY 

OBFLG 

OBMO 

OBYEAR 

OFLAG4 


p  and  angle  rejection  criteria 

M  TBINT 

e  - 

Flag  to  compute  mass  correction  analytically  or  by 
variant  ephemeris 

7905  (Prediction  Reliance) 

M 

Number  of  lines  to  be  output 
Pass  number  (Station  Prediction) 

Output  location  for (deg.) 

n  (BULGE) 
v 

vx 

V 

y 

V 

z 

Day  of  Observation 

Observation  flag  (What  quantities  were  observed) 
Month  of  observation 
Year  of  observation 
Residual  sets  to  output  flag 
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0FLAG5 

Set  up  by  OFLAG4  (Output  residuals  or  nol 

OLD  RMS  2 

Old  RMS  for  range  rate 

OLDRMS 

Old  rms  for  range  and  angles 

OLDUZ 

Old  U 

z 

OMEGAP 

Output  location  for  (jd  (deg.) 

0 

OPRTESQ 

1  +  J  1  -  e2 

ORMS2 

Old  rms  for  range  rate 

ORKS 

Old  rms  for  range  and  angles 

OUTAPE 

Output  tape  number 

PAGENO 

Pa.ge  number 

PAPRINT 

Output  location  for  Pao  (min.) 

PBUF 

P2  through  Pg 

(BULGE) 

PDAY 

Output  location  for  day 

PDAY 

Printing  location  for  day 

(IHEAD1) 

PFLAG 

Prediction  by  Station  pass  or  time  flag 

PHI 

0 

PISUN 

7T  =  4.929316613 

© 

PLSDE 

Plus  side  flag  (Station  Prediction) 

PN PRINT 

PPBUF 

P^  output  location 

P'  Buffer 

(BULGE) 

PPPFLAG 

Prediction  option  input 
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PP 

PREDBF 

PRTIME 

PRTIM 

PYEAR 

QQ 

RCNT2 

RCNT 

RDOTOR 

RDTFLG 

REJCNT 

REJFLG 

REV 

RGFLG 

RGSDL 

RHODT3 

RHOO 

RHO 

RMS  2 


+  1 

Buffer  to  save  new  epoch  time,  L,  a,  h 

Time  of  new  epoch 

Buffer  for  time  prediction 

Output  location  for  ye^r  (IHBAD1) 

Temporary  location 

Residual  count  for  p 

Residual  count  for  p  and  angles 


p  =0  flag  (Station  Prediction) 

Number  of  residuals  rejected 

Rejection  flag  if  any  observation  quantity  is 

rejected 

Revolution  number 
Range  residual  flag 
A p  (km) 

Storage  location  for  RHODT 

p  (8m/  3)  =  0.001225 

o  cm 

p  (DRAG) 

RMS  for  range  rate 
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RMSCHK 
RMS 
ROVA 
RPC ONI 
RPC ON 2 
RPCON3 

RPFLAG 

RPT 

RRFLG 

RRSDL 

RSFLAG 

RSFLG2 

RSFLG3 

SATEL 

SAIL 

SAVBUF 

SAVEO 

SAVE3 

SAVEM 

SENSNO 


50  km 

Root  mean  square  for  p,  and  angles 
r/a 

2e  =  0.03345100 
o 

5/4  e2  =  0.0003496775 

O 

„  „  .  Kg  radii 

F  =  y  P  A  — “ - r- 

o  'e  ,  .  2 

kemin 

Radiation  pressure  flag 
Number  of  passes  through  D.C. 

Range  rate  residual  flag 
A P>  (km/sec) 

Flag  for  CNTRL  for  intermediate  epoch 

First  time  through  CNTRL  for  intermediate  epoch 

Storage  location  WFLAG  (RESRET1) 

Satellite  number 

Output  location  for  Satellite  Number 

Buffer  for  saving  W  buffer  before  going  to  Runge- 
Kutta  back-up  routine 

Storage  location  for  cell  0 

Storage  location  for  cell  3 

^OBS’  -  V^OBS  '  £0BS 

Sensor  number  (Station  Prediction) 
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SETRD2 
SETRD 
SIGMA 1 

SIGMA 2 

SIGMA 3 

SIGMA 4 

SIGMA I 

SIGMA 

SINBOTH 

SINEPS 

SKNTRL 

SMLGR 

SOBFLG 

SOFI.AG5 

SV56 

SVDLTTO 

SVPRTIM 

SVRDTBD 

SVRHODT 

SWFLAG 


Address  after  last  cell  in  processed  observation  buffer 

Starting  address  of  processed  observation  buffer 
1 


5p"  for  P obs ’  ~p~p~  for  angles  observed 
a  (DRAG)  (TBINT) 

sin  (t jj+  e) 
sin  e  =  0.3978584 
Storage  location  for  KNTRL 
X  (DERIV) 

Storage  location  for  OBFLG 

Storage  location  for  0FLAG5 

Storage  location  for  index  registers  5  and  6 

Storage  location  for  DELTATO 

Storage  location  for  PRTIME 

Storage  location  for  RHODT 

Save  p  (Station  Prediction) 

Set  up  by  WFLAG 
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APPENDIX  C  (Continued) 


TCAP 

T  (Station  Prediction) 

TEMPI 

Temporary 

TEMP  2 

Temporary 

TEMP  3 

Temporary 

TEMP 

Storage  location  for  I'KTIME 

TERMS 

Least  Squares  Buffer 

TF 

tf  (final  output  time  -  time  pred) 

THETA 

e 

TLFLAG 

Point  of  closest  approach  flag  (Station  Predi 

impend 

Temporary  Ending  address  (Station  Prediction) 

TMPL0C1 

Temporary 

TMPLOC2 

Temporary 

TOPR 

Same  as  TO 

TOYPR 

Output  location  for  epoch  year 

T 

t  (minutes  since  epoch) 

TSI 

Beginning  time  point  (Time  Prediction) 

TS 

Ts  (°K) 

VFLAG) 

Variant  Ephemeris  Flag 
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APPENDIX  C  (Continued) 


WFLAG  Weight  Flag 

W  Fixed  or  Variable  integration  mode  (Runge-Kulta) 


WSUM2 

Normalized  RMS 

for  range  rate 

WSUM 

Normalized  RMS 

X10M0 

_1_ 

m0 

XBDGR 

. ' 

(BULGE) 

XDGR 

XB 

(DERIV) 

XDTGR 

•  \ 

XD 

(DRAG) 

XLAMD 

*E 

XLSNSBX 

L 

xo 

(RDPRES ) 

XLSNSBY 

L 

y® 

(RDPRES) 

XLSNSBZ 

L 

Z0 

(RDPRES) 

XLSUN 

L  (rad) 

a 

(RDPRES) 

XLSUNT 

t,  (rad) 

<9 

Km 

radii 

(RDPRES) 

XMPER 

XMSUN 

M  (rad) 

0 

(RDPRES) 

XNSUN1 

»  <d*S/H  >  • 

o  day 

0.9856473354  (THGRC) 

XNSUN 

n  (rad/  .  )  - 
0  mm 

0 .01 1 94638282  x  10'3  (RDPRES) 

XRDGR 

x'rp 

(RDPRES) 
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APPENDIX  C  (Continued) 


Y8DGR 

(BULGE) 

YDGR 

yv 

(DERIV ) 

YDTGR 

•  \ 
y  d 

(DRAG) 

YRDGR 

^RP 

(RDPRES ) 

ZBDGR 

.  \ 

ZB 

(BULGE) 

ZDGR 

.  \ 

2 

(DERIV) 

ZDTGR 

.  N 

ZD 

(DRAG) 

ZRDGR 

ZRP 

(RDPRES) 
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